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Turbulence and Skin Friction 


TH. von KARMAN, California Institute of Technology 
(Received November 3, 1933) 


1. HistorRICAL DATA 


HE first tests known on skin friction of flat 

surfaces in fluids are those carried out by 
M. G. Beaufoy! in 1793. About eighty years 
later W. Froude conducted a series of systematic 
experiments, investigating the influence of the 
velocity and the dimensions of the plate, espe- 
cially that of the length in the direction of the 
motion, on the skin friction. In 1883 O. Reynolds 
published his fundamental investigations on the 
similarity of fluid motion and introduced the 
dimensionless parameter known as Reynolds 
number (referred to hereinafter as R.N.). In 
1908 Lord Rayleigh pointed out that for geo- 
metrically similar arrangements, the skin friction 
coefficient, i.e., the friction per unit area divided 
by the head, must be a function of R.N. Suc- 
ceeding investigations of Zahm, Wieselsberger, 
Gibson, with air, of the British Froude National 
Tank, Gebers, Kempf and the Washington Navy 
Yard with water, were evaluated using the R.N. 
as parameter. 

Poiseuille gave in 1840 the theory of skin 
friction for the case of laminar flow in pipes. 
In 1850 G. Stokes published his calculations on 
the resistance of bodies in pure viscous flow, i.e., 
neglecting the inertia forces in the fluid. Prandtl 
and Blasius gave the mathematical theory of the 


' Beaufoy’s results were published by his son H. Beau- 
foy in 1834, 


laminar skin friction of plates in a free stream in 
1904 and 1907. The motion of the boundary 
layer introduced by Prandtl proved itself useful 
through the further developments of the theory 
both for laminar and turbulent flow. The first 
semi-empirical theory of turbulent skin friction 
was given in 1921 by Prandtl and by the present 
author. The recent development of the theory 
presented in this paper started with the author’s 
publications in 1930. 

The scientific and technical literature on skin 
friction is very elaborate. An excellent  bibli- 
ography, extending to 1930, has been compiled 
by A. F. Zahm and C. A. Ross.’ Additional 
references are given at the end of this paper. 


2. SKIN FRICTION AND MOMENTUM 


The modern theory of skin friction is based on 
the theorem that friction between a fluid and a 
solid is accompanied by an equivalent change of 
the momentum carried by the fluid. Let us 
consider the case of a fluid moving along a flat 
plate and compare the momentum carried by 
the fluid through a cross section upstream from 
the leading edge of the plate and through a 
cross section at the distance x downstream from 
the leading edge. We assume that the velocity 
is uniform over the first cross section and has 

2 Tentative bibliography on skin friction and boundary 
flow. Library of Congress, December, 1930. 
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the value U. The velocity normal to the second 
cross section at an arbitrary point therein shall 
be denoted by wu. The fluid mass passing in unit 
time through an element dS of the second cross 
section is evidently pudS, where p denotes the 
density of the fluid; this fluid mass had the 
momentum U per unit mass passing through 
the first cross section and has the momentum u 
per unit mass passing through the second cross 
section. Hence the loss of momentum of the fluid 
mass considered amounts to pu(U—u)dS and 
the total loss of momentum of the fluid in unit 
time to f pu(U—u)dS, the integral being taken 
over the cross section downstream. This amount 
is equal to the total frictional force acting on the 
portion of the plate extending from the leading 
edge to the distance x. | 

Let us consider a plate of the width } and 
introduce the coordinates x parallel to the 
direction of the main flow, y perpendicular to the 
flow and the plate with origin at the leading edge 


| 
U-u 
eu 
u TOTAL FORCE 


Fic. 1. Skin friction and boundary layer. 


(Fig. 1). Then the friction F, acting on the plate 
between x=0 (leading edge) and an arbitrary 
value of x will be 


pu(U—w)dy. (1) 


The friction acting on a strip bdx is obviously 
equal to dF; on the other hand, the same force 
can be expressed by the local friction per unit area 
or shearing stress at the wall 79 in the form drodx. 
Hence we obtain 


If (U (2) 
™=-—= pu(U—u)dy |. 


In all practical cases it is sufficiently exact to 
extend the integral from y=0 to a finite value 
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y = 6; we choose 6 in such a way that the velocity 
defect U—u is very small for values of y superior 
to 6. As a matter of fact, the influence of the 
friction is restricted to a narrow region along the 
plate; we call this region the ‘‘boundary layer”’ 
and 6 thickness of the boundary layer.” 
The choice of 6 is, of course, somewhat arbitrary, 
but as it will be seen, the results concerning the 
laws and the values of the skin friction are 
independent of the definition of 6. 

The consideration sketched above replaces the 
problem of the skin friction by the problem of 
the velocity distribution in the boundary layer. 
As a matter of fact, the friction itself can be 
determined as well by direct measurement as by 
measurement of the velocity distribution across 
the boundary layer. 


3. THE THEORY OF LAMINAR SKIN FRICTION 


In the case of laminar motion, the velocity 


distribution can be computed theoretically by. 


using the hydrodynamic equations in their some- 
what simplified form due to Prandtl (1904). The 
main result of this theory can be understood by 
using very little calculation if we assume similar 
velocity distributions over all cross sections of 
the boundary layer. Let us write, accordingly, 
u=U-f(y/5)=Uf(n) where n=y/6 is the “‘dis- 
tance ratio” referred to the thickness 6 of the 
boundary layer. The function f(y) equals zero 
for y=n=0 and is approximately equal to one 
for y=6 or n=1. The ‘momentum integral” 
pu(U—u)dy will be replaced by Sif 
-(1—f)dn. 

The integral /o'f-(1—f)dy is a pure number 
which shall be denoted by a. On the other hand, 
the friction at the wall 7, according to the 
general law of viscosity, equals the product of 
the viscosity coefficient » and the slope of 
velocity du/dy at the wall. Hence 


79 = w(du/dy)y—0= (3) 


The value (df/dn),~0 is also a pure number, we 
call it 8. Introducing the value of 7» and that of 
the momentum integral in Eq. (2) we obtain 


U d dé 
or 6—=--—:- 
6 dx d 


— 


A 


SKIN FRICTION 3 

: ‘ 152 — up distance from the leading edge and proportional 

By integration, 4 ” i Ua” to the geometrical mean value between viscosity 
and and density. 

5=[(28/a)(ux/pU)}}. (4) (c) The friction coefficients both for the local 


Putting this value of 6 in Eq. (3) we obtain the 
value of the local friction 


= uB(U/6) = (3a8)*- (on U8/x)! (5) 
and by integration, the value of the total friction 
F, = (6) 
We introduce the coefficient of the local friction 
(7) 

and that of the total (or mean) friction 
C;= F,/}+p: U?-b-x. (8) 


The mathematical theory gives for (2a8)} the 
numerical value 0.664. Hence with u/p=v 


c;=0.664/(Ux/v)} (9) 
and 
=1.328/(Ux/vr)}. (10) 


We call Ux/v the R.N. of the plate R,. 

Summarizing the results given by the equa- 
tions from (4) to (10), we obtain the following 
statements (cf. Fig. 2): 


Fic. 2. Distribution of laminar skin friction and thickness 
of laminar boundary layer. 


(a) The increase of the thickness of the lami- 
nar boundary layer 6 is proportional to the 
square root of the distance from the leading edge 
and inversely proportional to the square root of 
the relative velocity between plate and fluid. 

(b) The local friction per unit area is pro- 
portional to the 3/2 power of the velocity, 
inversely proportional to the square root of the 


and the mean friction are inversely proportional 
to the square root of the R.N. 


4. PowER LAws FOR TURBULENT 
SKIN FRICTION 


The method of the last section can be extended 
to the case of the turbulent boundary layer 
without investigating the actual mechanism of 
turbulent motion and by using merely the method 
of dimensional analysis. 

We assume, as before, that the velocity 
distribution is given by a function of the form 
u = Uf(y/6), but for the wall friction we use an 
empirical law of the form 1=}pU%g(U6/v). 
Using a ‘‘power law”’ as a convenient interpola- 
tion formula, we put 


= 3pU-const./(U6/v)™. (11) 


Evidently the momentum integral /o°pu(U—u)dy 
keeps the form apU*é, where @ is a numerical 
value, which depends on the shape of the velocity 
distribution in the boundary layer. Hence we 
obtain from the fundamental Eq. (2) 


apU?(dé/dx) =}pU*-const./(U6/v)™ (12) 
and 
6"(d6/dx) =const.-(v/U)™. (12a) 
Integrating (12a) and assuming 6=0 for x=0, 
we obtain 


5"+1/(m+1) =const.-(v/U)"-x 
or 
6=const.-x-(v/Ux)m/ (mt, (13) 


This equation determines the ‘development of 
the boundary layer” along the plate, if the 
“elementary friction law”’ (11) is known. 

Introducing the value of 6 in Eq. (11) for 70, 
we write 

or 
(14) 


We started with Eq. (11) which expressed the 
relation between the friction and the R.N. of the 
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Turbulence and Skin Friction 


TH. VON KARMAN, California Institute of Technology 
(Received November 3, 1933) 


1. HistorIcAL DATA 


HE first tests known on skin friction of flat 

surfaces in fluids are those carried out by 
M. G. Beaufoy' in 1793. About eighty years 
later W. Froude conducted a series of systematic 
experiments, investigating the influence of the 
velocity and the dimensions of the plate, espe- 
cially that of the length in the direction of the 
motion, on the skin friction. In 1883 O. Reynolds 
published his fundamental investigations on the 
similarity of fluid motion and introduced the 
dimensionless parameter known as_ Reynolds 
number (referred to hereinafter as R.N.). In 
1908 Lord Rayleigh pointed out that for geo- 
metrically similar arrangements, the skin friction 
coefficient, i.e., the friction per unit area divided 
by the head, must be a function of R.N. Suc- 
ceeding investigations of Zahm, Wieselsberger, 
Gibson, with air, of the British Froude National 
Tank, Gebers, Kempf and the Washington Navy 
Yard with water, were evaluated using the R.N. 
as parameter. 

Poiseuille gave in 1840 the theory of skin 
friction for the case of laminar flow in pipes. 
In 1850 G. Stokes published his calculations on 
the resistance of bodies in pure viscous flow, i.e., 
neglecting the inertia forces in the fluid. Prandtl 
and Blasius gave the mathematical theory of the 


' Beaufoy’s results were published by his son H. Beau- 
foy in 1834. 


laminar skin friction of plates in a free stream in 
1904 and 1907. The motion of the boundary 
layer introduced by Prandtl proved itself useful 
through the further developments of the theory 
both for laminar and turbulent flow. The first 
semi-empirical theory of turbulent skin friction 
was given in 1921 by Prandtl and by the present 
author. The recent development of the theory 
presented in this paper started with the author's 
publications in 1930. 

The scientific and technical literature on skin 
friction is very elaborate. An excellent  bibli- 
ography, extending to 1930, has been compiled 
by A. F. Zahm and C. A. Ross.? Additional 
references are given at the end of this paper. 


2. SKIN FRICTION AND MOMENTUM 


The modern theory of skin friction is based on 
the theorem that friction between a fluid and a 
solid is accompanied by an equivalent change of 
the momentum carried by the fluid. Let us 
consider the case of a fluid moving along a flat 
plate and compare the momentum carried by 
the fluid through a cross section upstream from 
the leading edge of the plate and through a 
cross section at the distance x downstream from 
the leading edge. We assume that the velocity 
is uniform over the first cross section and has 


2 Tentative bibliography on skin friction and boundary 
flow. Library of Congress, December, 1930. 
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the value U. The velocity normal to the second 
cross section at an arbitrary point therein shall 
be denoted by u. The fluid mass passing in unit 
time through an element dS of the second cross 
section is evidently pudS, where p denotes the 
density of the fluid; this fluid mass had the 
momentum U per unit mass passing through 
the first cross section and has the momentum u 
per unit mass passing through the second cross 
section. Hence the loss of momentum of the fluid 
mass considered amounts to pu(U—u)dS and 
the total loss of momentum of the fluid in unit 
time to f'pu(U—u)dS, the integral being taken 
over the cross section downstream. This amount 
is equal to the total frictional force acting on the 
portion of the plate extending from the leading 
edge to the distance x. 

Let us consider a plate of the width } and 
introduce the coordinates x parallel to the 
direction of the main flow, y perpendicular to the 
flow and the plate with origin at the leading edge 


U-u 


u TOTAL Force 


Fic. 1. Skin friction and boundary layer. 


(Fig. 1). Then the friction F, acting on the plate 
between x=0 (leading edge) and an arbitrary 
value of x will be 


pu(U—u)dy. (1) 
( 


) 


The friction acting on a strip bdx is obviously 
equal to dF; on the other hand, the same force 
can be expressed by the local friction per unit area 
or shearing stress at the wall 7» in the form brodx. 
Hence we obtain 


1dF d 
b dx dxt-/o 


In all practical cases it is sufficiently exact to 
extend the integral from y=0 to a finite value 
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y = 6; we choose 6 in such a way that the velocity 
defect U—u is very small for values of y superior 
to 6. As a matter of fact, the influence of the 
friction is restricted to a narrow region along the 
plate; we call this region the “boundary layer”’ 
and 6 thickness of the boundary layer.” 
The choice of 6 is, of course, somewhat arbitrary, 
but as it will be seen, the results concerning the 
laws and the values of the skin friction are 
independent of the definition of 6. 

The consideration sketched above replaces the 
problem of the skin friction by the problem of 
the velocity distribution in the boundary layer. 
As a matter of fact, the friction itself can be 
determined as well by direct measurement as by 
measurement of the velocity distribution across 
the boundary layer. 


3. THe THEORY OF LAMINAR SKIN FRICTION 


In the case of laminar motion, the velocity 
distribution can be computed theoretically by 
using the hydrodynamic equations in their some- 
what simplified form due to Prandtl (1904). The 
main result of this theory can be understood by 
using very little calculation if we assume similar 
velocity distributions over all cross sections of 
the boundary layer. Let us write, accordingly, 
u=U-f(y/6)=Uf(n) where n=y/6 is the “‘dis- 
tance ratio” referred to the thickness 6 of the 
boundary layer. The function f(y) equals zero 
for y=n=0 and is approximately equal to one 


-for y=6 or n=1. The ‘momentum integral” 


So’ pu(U—u)dy will be replaced by pU6 Sif 
-(1—f)dn. 

The integral fi'f-(1—f)dn is a pure number 
which shall be denoted by a. On the other hand, 
the friction at the wall 7, according to the 
general law of viscosity, equals the product of 
the viscosity coefficient » and thes slope of 
velocity du/dy at the wall. Hence 


ro = 


The value (df/dn),~0 is also a pure number, we 
call it 8. Introducing the value of 7) and that of 
the momentum integral in Eq. (2) we obtain 


U d dé 
or 
6 dx dx pUea 
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up 
By integration, 36? =———-x 
pUa 
and 
°6=[(28/a)(ux/pU) ]}. (4) 


Putting this value of 6 in Eq. (3) we obtain the 
value of the local friction 


ro= uB(U/0) (5) 
and by integration, the value of the total friction 
F, = (2a8)*+ (puxU*) (6) 
We introduce the coefficient of the local friction 
Cs = 19/3 pU? (7) 

and that of the toval (or mean) friction 
C,= F,/}-p- U?-b-x. (8) 


The mathematical theory gives for (2a)! the 
numerical value 0.664. Hence with uw/p=v 


=0.664/(Ux/v)! (9) 
and 
Cy =1.328/(Ux/v)}. (10) 


We call Ux/v the R.N. of the plate R,. 

Summarizing the results given by the equa- 
tions from (4) to (10), we obtain the following 
statements (cf. Fig. 2): 


Fic. 2. Distribution of laminar skin friction and) thickness 
of laminar boundary layer. 


(a) The increase of the thickness of the lami- 
nar boundary layer 6 is proportional to the 
square root of the distance from the leading edge 


. and inversely proportional to the square root of 


the relative velocity between plate and fluid. 
(b) The local friction per unit area is pro- 

portional to the 3/2 power of the velocity, 

inversely proportional to the square root of the 


distance from the leading edge and proportional 
to the geometrical mean value between viscosity 
and density. 

(c) The friction coefficients both for the local 
and the mean friction are inversely proportional 
to the square root of the R.N. 


4. PowER LAWS FOR TURBULENT 
SKIN FRICTION 


The method of the last section can be extended 
to the case of the turbulent boundary layer 
without investigating the actual mechanism of 
turbulent motion and by using merely the method 
of dimensional analysis. 

We assume, as before, that the velocity 
distribution is given by a function of the form 
u = Uf(y/6), but for the wall friction we use an 
empirical law of the form 1t)=}3pU*g(U6/r). 
Using a ‘‘power law” as a convenient interpola- 
tion formula, we put 


t) = $pU?-const./(U6/v)™. (11) 


Evidently the momentum integral /j°pu(U—u)dy 
keeps the form apU*5, where a is a numerical 
value, which depends on the shape of the velocity 
distribution in the boundary layer. Hence we 
obtain from the fundamental Eq. (2) 


ap U?(dé/dx) =}pU*-const./(U6/v)™ (12) 
and 
6"(d6/dx) =const.-(v/U)™. (12a) 
Integrating (12a) and assuming 6=0 for x=0, 
we obtain 
6"+!/(m+1) =const.-(v/ U)"-x 
or 


6=const.-x-(v/ (13) 


This equation determines the “development of 


gine boundary layer’ along the plate, if the 


‘elementary friction law”’ (11) is known. 
Introducing the value of 6 in Eq. (11) for 7o, 
we write 


or 
(14) 


We started with Eq. (11) which expressed the 
relation between the friction and the R.N. of the 
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boundary layer. It is logical to assume that this 
relation is similar to that obtained for a channel 
or pipe with constant cross section. For laminar 
flow we have therefore, according to Poiseuille’s 
law, m=1 and therefore m/m+1=4 so that Eq. 
(13) leads to the results of the last section. For 
turbulent flow tests on pressure drop in smooth 
pipes give in the range 3000< Ur/v<70,000 
(ry =radius of circular pipe) good accordance with 
the empirical Eq. (11) choosing the exponent 
m=1, 
present author obtained from (14) the formula 
for the turbulent skin friction coefficient of a 


smooth plate: 
(15) 


and with the value of the numerical constant 
deduced from the pipe tests 


It is easy to show that the coefficient of the 
mean friction is equal to C;=c;-(m+1)/m and 
therefore 


C; =0.074/(Ux/v)"*. (16) 


This formula checks the tests with smooth plates 
in the range extending to about Ux /v=3-10°. 
Extension of the tests to larger R.N. made it 
evident that the power formula with a constant 
exponent could not serve for a larger range of 
R.N., either in the case of the pipe or in the case 
of the plate. It was found that the exponent m 
decreases with increasing R.N. This proved the 
assumption of similar velocity distributions to 
be incorrect. In order to find a law valid in the 
whole range and especially valid for large R.N. 
allowing extrapolation to full scale conditions, a 
more systematic investigation on the nature of 
turbulent flow became inevitable. This new 
development of the theory achieved in the last 
few years is shortly reviewed in the following 
sections. 


5. Tue MECHANISM OF TURBULENT FLOW 


In laminar flow the‘interaction between ad- 
jacent fluid layers consists of so-called ‘‘molecular 
friction’; the shearing stress, i.e., the friction 
per unit area, is equal to the product of the rate 
of shear (slope of velocity in the case of parallel 
flow) and the viscosity coefficient. The latter is a 


Using this exponent, Prandtl and the 
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characteristic physical parameter of the fluid, 
dependent chiefly on the temperature, practically 
independent of the velocity. In turbulent flow 
an additional interaction is due to the momentum 
transfer from layer to layer due to the velocity 
fluctuations. Let us consider a flow, consisting 
of parallel mean motion in the x-direction, the 
temporal mean value u of the velocity being a 
function of y only. O. Reynolds has shown that 
the momentum transport in unit time and 
through unit area due to the fluctuations repre- 
sents an apparent shearing stress for the mean 
motion. This apparent shearing stress, or “‘turbu- 
lent friction,” has the magnitude += —pu’v’ 
where w’ and v’ are the fluctuations of the velocity 
components in the x- and y-direction, and the 
dash denotes temporal mean value of the product. 
In fact the fluid volume passing through a unit 
area perpendicular to the y-direction is equal to 
v’ and the excess of momentum in the direction 
of the mean flow carried by the fluid per unit 
volume is equal to pu’. The negative sign 
corresponds to the definition of positive friction 
in such a way that 7 is considered as positive if 
the fluid layer at the distance y from the wall is 
accelerated by the outside flow; accordingly, the 
momentum transfer produces friction of positive 
amount if the fluid portions moving toward the 
wall (i.e., having negative v’ components) carry 
positive excess of momentum. 

In the case of velocity fluctuations with a 
random distribution of the «’ and v’ components, 
the product u’v’ has the mean value zero. The 
shearing stress is different from zero only in the 
case of a certain correlation between wu’ and of 
It is easy to see that a correlation can be expected 
in all cases in which the mean velocity w is 
variable with y. Let us assume for instance that 
u(y) is increasing with y (Fig. 3). Then the 
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Fic. 3. Turbulent exchange of momentum. 
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fluid masses moving in the negative y-direction 
(v'<0), ie., coming from a region of higher 
velocity, probably have larger u-components 
than the masses with positive v’ values, i.e., 
coming from a region of lower velocity. Now 
evidently larger wu-velocities than the mean 
velocity at the place considered appear as posi- 
tive, smaller u-velocities as negative w’ fluctua- 
tions; hence negative v’ values will be connected 
with positive wu’ values and vice versa. Therefore, 
we have to expect a negative value for the mean 
product «’v’ and a positive value for the shearing 
stress. Fig. 4 shows schematically the random 
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Fic. 4. Diagram of velocity fluctuations. 


and correlated distribution of the velocity fluctu- 
ations. 

The main characteristics of the turbulent flow 
at a certain point are the magnitude of the 
fluctuations wu’, v’, w’, and the correlations be- 
tween them. The theory of correlations uses a 
certain ratio between mean values as ‘correlation 
factor.’’ For instance, the correlation between 
the components _«’ and v’ is measured by the 
ratio where (uw)! and (v”)! are 
the square roots of the mean squares of the 
particular components. This correlation factor 
is zero for random distribution and is equal to 
unity if uw’ and v’ are proportional. The ellipse 
shown in Fig. 4 becomes a circle in the first case 
and a straight line in the latter. 

Considering the ‘‘flow pattern” in the neigh- 
borhood of a certain point, we have to introduce 
a further characteristic parameter of the turbu- 
lent flow: namely a length which is characteristic 
for the size of the region involved in the turbulent 
exchange, i.e., in the turbulent momentum trans- 
fer. The corresponding parameter in the molecu- 
lar mechanism of friction is the mean free path. 
Prandtl uses for the length-parameter of the 
turbulent friction-mechanism the expression 
“mixing length.” His particular picture of the 


turbulent exchange is based on the assumption 
that the fluid masses displaced perpendicular to 
the direction of mean flow carry their momentum 
invariably over a certain length perpendicular 
to the mean flow. The actual mechanism may 
be more complicated; in any case the conception 
of a characteristic length is necessary for the 
understanding of the fundamental laws of turbu- 
lent flow. Let us compare, for instance, the 
turbulent fluctuations in the boundary layer of 
an airfoil placed in the wind tunnel and those in 
the boundary layer produced by the natural 
wind on the ground. Mean velocities and ve- 
locity fluctuations may have the same order of 
magnitude; it is obvious that the length scale of 
the flow pattern in the two cases will be very 
different. 

There is a further factor to be considered in 
connection with the characteristic length or with 
the length scale of the flow pattern, namely the 
time-factor or the time-scale of the fluctuations. 
“Large scale turbulence” will be connected in 
general with fluctuations of large period, and 
“small scale turbulence’”’ with fast fluctuations. 

The characteristic features of turbulent flow 
have been dealt with in some detail and from a 
general point of view, because the author feels 
that clear notions on this question are of im- 
portance for understanding the phenomena and 
for further research. Common language, even in 
technical papers, uses the term “turbulence” in a 
rather indefinite way, including in it for instance 
“vortex motion” in general. The scientific term 
“turbulence,”” as used in this paper, implies 
“irregular fluctuations,’’ governed by laws of 
some statistical equilibrium. To be sure, turbu- 
lence may be originated by vortices as in the 
case represented in Fig. 5, showing the flow 
through a honeycomb; but we call such a 
motion “‘turbulent” only if the regular pattern, 
because of the intermingling of a great number 
of vortices, disappears, as happens farther down- 
stream from the honeycomb. Similarly, vortices 
produced by obstacles at the ground contribute 
essentially to atmospheric turbulence, but also 
in this case we restrict the term turbulence to the 
statistical phenomenon of the mass exchange and 
keep this case apart from other conceptions, such 
as individual vortices (tornados, vortex sheets, 
etc.) or regular atmospheric waves. 
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Fic. 5. Vortex rows produced by a 


6. THE Basic LAWs FOR THE MEAN VELOCITY 
DISTRIBUTION IN TURBULENT FLOW .- 


Adding the momentum transfer by turbulent 
fluctuations to the laminar friction, O. Reynolds 
obtained the following expression for the total 
friction (shearing stress) : 


r=y-du/dy—pu'v’. (17) 
Calculating + from the measured values of the 
pressure drop in pipes or the skin friction of 
plates, the values yu(du/dy) from the measured 
mean velocity distribution, we find that in the 
practically important range of high R.N. the 
laminar friction is a negligible part of the total 
friction except in the immediate neighborhood of 
the wall. Hence we conclude that except for 
this region the turbulent exchange is practically 
independent of the viscosity of the fluid. This 
assumption is essentially substantiated by experi- 
mental facts concerning the velocity distribution 
in pipes and channels. Let us consider, for 
instance, turbulent flow between two parallel 
walls (rectangular pipe with great aspect ratio). 
The velocity in the center will be denoted by U, 


honeycomb. (Photograph by Ahlborn.) 
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the velocity at a distance y from the wall by u, 
the distance between center and wall by 7, the 
wall friction by ro. The difference (U—u) will 
be called the “‘velocity defect.’’ We find experi- 
mentally that for high R.N. the ratio (U—u) 
(ro/p)! is a general function of the ratio y/r or 
(r—y)/r, i.e., independent of R.N. and also 
independent of the roughness of the wall. 

A simple interpretation of this result can be 
given by comparing cases in which both 79 and r 
have the same values, but the roughness of the 
wall is different (Fig. 6). The velocity defect 
(U—u) is the relative velocity between the fluid 
moving in the center layer and the fluid moving 
in a layer at a distanc€é (r—y) from the center. 


y Rough Smooth Wal/ 
~U-u 
r 
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Channel Center | 


Fic. 6. Comparison of velocity distributions in a smooth 
and in a rough channel. 
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Considering the flow relative to the fluid moving 
in the center layer, the fact that (U—n) is 
independent of the roughness of the wall and is 
proportional to (79/p) means that the action of 
the wall friction is independent of the way the 
friction is produced: whether by slow displace- 
ment of a very rough wall or by fast displacement 
of a smooth wall. The mechanism of the turbu- 
lent exchange seems to be the same in both 
cases, except in the immediate neighborhood of 
the walls. 

The same experimental result is true for the 
flow in a pipe with circular cross section. De- 
noting the radius of the pipe by 7 we obtain 


(U—u)/(70 p)i=fly r), (18) 


where f(y/r) is a universal function. 

The quantity (7o/p)? has evidently the di- 
mension of a velocity. In aerodynamics it is usual 
to characterize a velocity by the corresponding 
pressure head; in the inverse way we can charac- 
terize a certain value of the pressure p by the 
corresponding velocity (2p/p)?. In a_ similar 
way we can assign the velocity (279/p)' to the 
shearing stress ro. The pressure head corre- 
sponding to this velocity is numerically equal to 
the shearing stress. Instead of (279/p)! we will 
use (to/p)*. To be sure it would be more logical 
to use the former quantity, but (79/p)! has been 
used already in so many publications that a 
change would likely produce confusion. We call 
the quantity (7o/p)* the friction velocity. The 
ratio y/r can be called the relative distance from 
the center. 

Hence the statement expressed by Eq. (18) 
can be worded in the following way: The ratio 
( velocity defect 

friction velocity 
relative distance from the center. 

Let us now consider the velocity distribution 
relative to the wall. At the wall itself the velocity 
is zero or equal to the velocity of the wall, if the 
wall is moving. Hence, at the wall the velocity 
components w’ and v’ vanish, and the momentum 
transfer is zero, so that the friction is merely 
laminar friction. Thus we obtain a “laminar 
sublayer’’ adjacent o the wall with gradual 
transition to the turb ilent region. Therefore, the 
velocity distribution 1s function of the distance 


) is a universal function of the 
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y will contain not only the wall friction 7> and 
the density p, but also the kinematic viscosity v. 
Let us restrict the investigation to the case of 
perfectly smooth walls. In this case it can be 
expected that near the wall wu is in fact fully 
determined by to, y, p and v. Dimensional 
analysis shows easily that the only possible 
combination is given by the formula 


u=(To p)*g((to/p)*y ‘v), (19) 


where g is another universal function. 

The parameter (70/p)*y v, introduced first by 
L. Prandtl, is built up in the same way as the 
R.N.; it contains the friction velocity, the 
distance and the kinematic viscosity. We call it 
the friction-distance parameter. 
velocity 


Hence Eq. (19) says the ratio —— : 
friction velocity 


near the wall is a universal function of the friction- 
distance parameter. 


7. SIMILARITY OF THE TURBULENT 
FLow PATTERN 


In order to determine the function f(y 7) 
representing the velocity distribution over the 
cross section of a pipe of a channel, the author 
introduced the assumption that the flow pattern 
of the turbulent interchange in the neighborhood 
of different points is similar and only different 
as far as length and time scale are concerned. 
This assumption appears to be justified in all 
cases in which the characteristic length of the 
turbulent exchange is small in comparison with 
the dimensions of the cross section. 

The following conclusions are deduced from 
the above-mentioned assumption: 

(a) The velocity fluctuations are proportional 
to a characteristic length / of the flow pattern 
and to the slope of the mean velocity du/dy. 

(b) The shearing stress 7, variable over the 
cross section, is at every point proportional to 
pl?(du/dy)*. Hence / is heretofore only a relative 
measure of the flow pattern; we can absorb the 
proportionality factor and define / by writing 
=pl*(du/dy)?. 

(c) The characteristic length is proportional 
to the ratio |(du/dy)/(d?u/dy*)|. We write 
(du,dy)/(d’u/dy*)|, where « is a numerical 
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Fic. 7. Correlation between turbulent velocity fluctuations. 
(a) Wattendorf (C.1.T.). (b) Reichardt (Géttingen). 


factor. It will be seen that the constant « occurs 
in most of the formulae given in the later 
sections; it appears as a “‘universal constant”’ of 
the turbulent exchange. 

The formula 7 = pl?(du/dy)? was first given by 
Prandtl in his theory of the mixing length, 
assuming a special kind of mechanism for the 
turbulent exchange of momentum. 

The assumption of similar flow patterns is 
identical with the assumption of constant corre- 
lation between the components of the velocity 
fluctuations. There is not much experimental 
evidence on the correlation of the fluctuations, 
but in some cases the magnitude of the fluctua- 
tions has been measured. Recently F. Wattendorf 
of the California Institute of Technology meas- 
ured the value of wu” over the cross section of a 
rectangular channel with aspect ratio 1 : 18, 
using a hot wire anemometer suitable for fast 
fluctuations. In this case the value of u’v’ could 
be deduced from the measurement of the pressure 
drop, assuming += —pu’v’ and taking into con- 
sideration that the shearing stress 7 in such a 
channel is almost exactly proportional to the 
distance from the center. Thus the ratio u’v’/u” 
could be calculated; the results are given in Fig. 
7 (curve (a)). The points connected by the curve 
(6) correspond to similar measurements by 
Reichardt (Géttingen). The curves show the 
following results: ts 

(a) The ratio w’v’/u” is almost constant over 
the main part of the cross section, corresponding 
to constant correlation. 

(b) For the center w’v’/u’” 0; this is evident 
by symmetry considerations, since in the center 
no reason for correlation exists. Accordingly, in 
the center part the similarity assumption cannot 
be correct. 
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(c) The curve (6) drops near the wall. How- 
ever, Reichardt’s experiments are carried out 
at a relatively low R.N., so that the range of 
viscosity influence is comparatively large. The 
estimated range in which pwdu/dy cannot be 
neglected in comparison with —pu’v’ is indicated 
by the dotted lines for both cases. The way we 
can estimate the thickness of the range influenced 
by laminar friction is shown in the latter part 
of this section. . 

The assumption of similarity of the turbulent 
flow pattern gives the following formula for the 
velocity distribution: 


U-u 1 y\! 
(To p)? K Yr 


The formula checks the measured velocity 
distributions very well (Fig. 8). 
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Fic. 8. Measured velocity dis:ributions compared with 
the theoretic curve. 
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Fic. 9. Universal velocity distribution near a smooth wall. 


For small values of y/r, we obtain 


Umax — tl —(1/k)(70 ‘p)} log (y 


or (20) 


log yl 


This result can be seen directly in the following 
way: 

Let us consider the idealized case that a 
constant shearing stress is transferred in a 
parallel flow along the wall. In this case an 
exact similarity of the flow pattern can be 
expected; in fact looking at the flow from the 
point of view of two observers travelling with 
the fluid with a velocity equal to the mean 
velocity of the fluid, at two different distances 
vy, and ye from the wall, nothing can be different 
in the flow pictures, conceived by the two 
observers, except the length and scale of the 
flow pattern. Thus the characteristic length of 
the flow pattern is proportional to y. The 
distribution of the mean velocity is, according to 
Taylor’s theorem, given by 


u(y+ Ay) =u(y)+u'Ayt+ iu" Ay?+---. 


It is evident that the similarity of the flow 


pattern involves =u - There- 


fore we put 


» 2 9 


and obtain 
u log y). (21) 


Comparing Eq. (21) with the general Eq. (19), 
we obtain 


u=(rto/p)*La+(1/x) log p)'y/v) (22) 


This equation was given by the present author 
in 1930. The numerical value found by the author 
for the universal constant is 0.38. Nikuradse 
found later that x=0.40 fit the experimental 
values better. In Fig. 9 are plotted Nikuradse’s 
velocity distribution measurements — using 
u+=u/(to/p)' and logio yt =logio ((t0/p)*y/r) 
as coordinates. For y*<110 individual points 
are shown, for higher values of y* a certain 
number of measured points are replaced by their 
centers of gravity. It is seen that the equation 
ut =5.5+5.75 logiy y* fits the measurements very 
well. The factor 5.75 corresponds to «=0.40. 
There is a systematic deviation for values of 
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yt =(10/p)'y/v<30. Obviously for small values 
of the friction distance the similarity assumption 
leading to the logarithmic law cannot be correct. 
We remember that at the wall itself the flow is 
laminar. If we assume a pure laminar sublayer 
corresponding to the constant value 7, of the 
friction, the corresponding velocity distribution 
would be given by “=(r0/mu)-y or 


ut=yt. (23) 


The two lines representing (22) and (23) meet 
at the point (70/p)'y/v=11.5. Hence if there 
were an abrupt transition between pure laminar 
and turbulent regions, the ‘‘thickness of the 
laminar sublayer’’ would be given by y*=11.5 
or y=6,;=11.5v/(70/p)'. Actually the transition 
is gradual and the value y*=30 appears as a 
reasonable limit for the region influenced by 
laminar friction. 

Before the author found Eq. (22), power laws 
were used for Eq. (19). If we put 


(24) 


where A is a numerical constant, it is easy to see 
that a relation must prevail between the expo- 
nent m of the velocity distribution and the 
exponent m in the friction law, mentioned in 
Section 4. In fact if we solve Eq. (24) for ro, 
we obtain 

1 


(25) 


To=dpu?- : 
(u y / 


Therefore m=2n/(n+1), or n=m/(2—m). For 
m=1/4 we obtain n=1/7. We mentioned that 
for large R.N. m diminishes, as also does n. In 
Fig. 9 is shown how far the power law formulae 
with n»=1/7 and 1/8 represent an approximation 
to the logarithmic formula in certain ranges. 


8. SKIN FRICTION IN SMOOTH PIPES 


Combining Eqs. (18) and (22), i.e., the de- 
velopments for the velocity distribution starting 
from the center and from the wall, the author 
obtained in 1930 a rather simple formula con- 
necting the frictional loss in pipes and R.N., the 
jirst formula which at the same time is based on 
theory and checks the experiments. 

Let us consider the velocity u at a point y 
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outside of the laminar range, but near enough 
to the wall so that the logarithmic formula (22) 
for the velocity distribution is sufficiently exact. 
In such a region the two Eqs. (18) and (22) 
overlap, so that both are valid. Hence we write 


u 1 (ro/p)*y 
- =a+— log, —, (26) 
(to/p)? K 
and at the same time 
U-u r—y 
-1(—). (27 
(to/p)* 


If both equations hold, Eq. (27) for small values 
of y, must have the form: 
U-u 1 
-= const. —— log -, 
(t0/p)! KF 


(28) 


and adding (26) and (28), we obtain 


U 1 
—=const.+-— log (29) 


(to/p)! K 


Let us introduce the friction coefficient referred 
to the maximum velocity by the relation 


= pU? (30) 
and the R.N. referred to the same velocity 
R=Ur/v. (31) 
Then from Eq. (30) follows immediately 
(2/c;)}=const.+ (1/x) log. (32) 


Eq. (32) can easily be checked by plotting the 
values of 1/c;} as function of log (R(c;)'). If the 
equation holds, 1/c;} appears as a linear function 
of log (R(c,)'), and the slope of the straight line 
representing 1/c;' as a function of log (R(c;)*) 
determines the value of the universal constant x. 
This plotting is shown in Fig. 10. The circles 
represent centers of gravity of groups of Niku- 
radse’s measured points. The crosses are repre- 
sentative of mean values taken from measure- 
ments of Stanton and Pannell. The straight line 
corresponds to the equation 


1/c;§=3.60+4.15 log io (R(c,) (33) 


The factor 4.15 in the last term corresponds to 
the value of the constant «=0.39. 
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Fic 10. Plotting of 1/cy? against logio (R(c;)*) for 
smooth pipes. 


It is easy to show that a formula of the same 
form with another additive constant is valid for 
the skin friction coefficient referred to the mean 
velocity U instead of to the maximum velocity 
U. This remark is due to L. Prandtl. 

According to Eq. (27), with the exception of a 
small region near to the wall, (U—1)/(r»/p)! is 
a universal function of y/r, independent of the 
R.N.; integrating (U—w) over the cross section 
we obtain 

(U-—u)dA =(U-—U)-A 
so that 
(U—U),/(ro/p)! =const. 
or 
1 (34) 


On the other hand the value of R(c;)) is 
invariable if both R and c, are referred to the 
same parameter, i.e., both to the maximum or 
both to the mean velocity. 

The friction coefficient \ used in the hydraulic 
practice is defined by the relation 


h=(1/d)(U?/2¢), 


where h denotes the pressure loss measured in 
height, / the length, d the diameter of the pipe. 
The factor \ is the double of the coefficient Cy. 
Nikuradse and Prandt! gave, following the 
author’s deductions, the formula 


= —0.84+2 login 


The factor 2 in this formula corresponds to the 
value 0.403 for the universal constant x. The 
value of « can be determined from (a) velocity 
distribution measurements, (b) tests in smooth 
pipes, (c) tests in rough pipes, (d) friction 
measurements on plates. Nikuradse gives formu- 
lae identical in form to those given by the 
author previously for the three cases (a)—(c), 
but he determines the numerical coefficients 
independently for each case. Comparing his 
formulae, we find «=0.40 for (a), «=0.403 for 
(b) and (c). Plotting his results and those of 
Stanton and Pannell and Fromm, we find that 
x=0.39 fits both items (b) and (c) well; the 
same constant is right for the skin friction of 
plates (d). For the velocity distribution near to 
the wall «=0.40 checks the experiments better, 
for the velocity distribution in greater distance 
from the wall «=0.38 fits better. In view of the 
scattering of the measured values and the quite 
considerable differences between the measure- 
ments of different observers, it seems to the 
author that the accordance between the x values 
obtained in different ways is quite satisfactory 
and substantiates the existence of such a uni- 
versal constant entering in all four problems 
mentioned. 


9. SKIN FRICTION OF SMOOTH PLATES 


The results obtained and checked by com- 
parison with pressure drop measurements have 
been applied by the author to the computation 
of the skin friction of smooth plates. The basic 
assumption is that Eq. (32) holds for the turbu- 
lent boundary layer, if the R.N. involved in the 
equation is replaced by the R.N. of the boundary 
layer, Rs=U6/v. The calculation leads to a 
formula of the form 


2! /c;}=const.+(1/x) log (R,c,) (35) 
both for the coefficient of the local friction 
70/3 pU? (36) 
and for the coefficient of the total or mean 
friction 
(37) 


where A is the surface area of the portion 
considered. The R.N. is referred to the velocity 
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U outside of the boundary layer and the length 
x of the portion of the plate considered. 

The Eq. (35) settles some questions often 
discussed without decisive answer in previous 
technical papers on skin friction. It shows for 
instance that the contribution of equal areas to 
the skin friction of a plate diminishes slowly but 
indefinitely with increasing distance from the 
leading edge. Also the coefficient of the mean 
friction diminishes with increasing R.N., i.e., 
with increasing velocity and increasing length. 
Formerly it was assumed that both the local 
and the mean friction coefficients have certain 
finite inferior limits. 

Concerning the numerical values of the con- 
stants involved in Eq. (35), the same method of 
plotting can be used as shown above for the 
case of the pipes. In Fig. 11 are plotted the 


30; 
25 
20 
15 
105 3 4 5 6 7 


log,(R 


Fic. 11. Plotting of 1/cy} against logio (Rey) for 
smooth plates. 


values of the local friction coefficients according 
to measurements of G. Kempf. These experi- 
ments are very remarkable because the values 
of the local friction were obtained by direct 
measurements of the force acting on small 
movable plates arranged at different points of a 
very long pontoon. 

The equation of the straight line representing 
the measurements is 


logio (R-C)). (38) 


E. Schoenherr found that the coefficient of the 
mean friction C; is given by 


0.242 C;} = logio (R.C;), (39) 


KARMAN 


corresponding to x =0.395, or 
1 Cy? =4.15 logio 


Between the coefficients of the local and the 
mean friction, Schoenherr found the relation 


c, =0.558C,/0.558+2(C,)}, (40) 


where the numerical constant 0.558 represents 
«(2)}. 

The calculation leading to Eq. (35) can be 
found in the author’s paper presented at the 
International Congress for Applied Mechanics 
(Stockholm 1930) or in his paper on the theory 
of skin friction in the Proceedings of the Ham- 
burg conference on problems of ship propulsion 
(1932). The Eq. (35) is obtained by neglecting 
terms of the order C;) in comparison with 
1/C,;*. For moderate values of the R.N. the 
terms omitted are not quite negligible, but it 
appears that by slight changes in the constants 
of the Eq. (35) the whole practical range can be 
covered with sufficient exactitude by this simpler 
formula. Prandtl repeated the calculation in a 
slightly different way and gave the results of 
his calculation in the Gottinger Ergebnisse of 
1932. His results are very similar to those ob- 
tained by the present author. 

Eqs. (38) and (39) have the disadvantage that 
they cannot be analytically resolved for c; and 
C;. The Table I gives the values of cy; and C; as 
function of R. 


TABLE I. 

R, Cy R, Cy 
2x 108 0.00614 0.00480 0.00263 0.00222 
3 562 444 | 3 247 216 
4 529 419 | 4 237 210 
5 506 402 | 5 229 196 
7 473 379 | 7 218 187 

10° 441 357 |108 207 179 
1.5 108 408 332 | 2108 189 163 
2 387 317 3 179 156 
3 360 296 | 5 167 145 
4 342 283 |10° 155 135 
5 330 274 | 2x10° 141 125 
7 att 200 5 127 113 

107 293 245 |10'9 118 105 


Schlichting suggested that the Eq. (39) can 
be replaced by an empirical formula C;=0.455 
(logio R,)?**. This formula replaces the Eq. (39) 
fairly well in the range 10°<R,< 10°. 
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Fic. 12. Skin friction of smooth plates compared with the theory. 


10. DiscUssION OF THE EXPERIMENTs. TRAN- 
SITION BETWEEN LAMINAR AND 
TURBULENT REGION 
Fig. 12 represents the results of the most 
important test on skin friction of smooth plates. 
The experimental points are, in general, replaced 
by curves connecting the mean values. In 
addition, the theoretical curve for the coefficient 
of the laminar skin friction and the curve for 
turbulent skin friction corresponding to Eq. (39) 

are drawn. 

The tests of Wieselsberger, Zahm and Gibson 
were wind tunnel tests on stationary models; all 
other tests were made by towing flat plates in 
water. In all tests except those made by G. 
Kempf the total force acting on the model plate 
was measured and a correction was made for the 
additional drag due to the finite thickness of the 
plate. In Kempf’s tests the local friction was 
measured and the total friction computed by 
integration. 

It is interesting that in the range 2X10°<R 
<10° the results of the wind tunnel tests check 
with the theoretical values for turbulent friction 
and the tests of the Froude national tank check 
approximately with the theoretical values for 
laminar friction. Geber’s results in the range 
10°<R<10' show transitions between the two 


theoretical curves. 


The transition between laminar and turbulent 
boundary layers along a flat plate was investi- 
gated by Burgers, Hegge Zynen, Hensen, Fage, 
and recently by Schoenherr. The ideal case of a 
mathematically thin plate can be approximated 
by using plates with sharp leading edges. It is 
to be expected that under similar conditions the 
transition depends on the R.N. of the boundary 
layer R; = Ué/v. Hegge Zynen and Burgers found 
that in their wind tunnel experiments the critical 
value of R; was between 1650 and 3500; Hansen 
found the value R;, = 3100. It appears that there 
is no definite characteristic critical value of R; 
for the transition; the value R;,. depends on the 
amount of perturbations, for instance on the 
turbulence in the outside stream. Schoenherr 
found that the transition was substantially 
facilitated by producing artificial turbulence in 
the water upstream from the leading edge. He 
obtained the transition curves (b) and (c) shown 
in Fig. 13 under normal conditions and the 
curve (a) for artificial turbulence produced by 
roughening the forward edge of the plates. He 
also found that if artificial turbulence is produced 
in the water upstream from the leading edge the 
skin friction values, from R=2 10° on, fall on 
the curve for pure turbulent friction, as in the 
tests of Zahm and Wieselsberger. In general, 
rounding the leading edge has a similar influence. 
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Fic. 13. Transition between laminar an 


The range of critical values of R; can be 
assumed to be between 1600 and 6000. 

Assuming that the transition takes place 
abruptly at a certain value of R; a family of 
“transition curves” for C; can be calculated as 
function of the R.M. of the plate with the 
critical value R;, as parameter. Such curves are 
drawn in Fig. 13. Curves I and II correspond 
to somewhat different assumptions. In both cases 
a laminar layer is assumed, until a certain value 
of R; is reached. The turbulent layer beyond the 
transition point is calculated in case I, as if it 
had started with the thickness zero at the 
leading edge; in case II it is assumed that at the 
transition point the momentum carried by the 
boundary layer remains unchanged. Geber’s 
experimental curves fit the calculated curves 
quite well. However, in view of the Schoenherr 
curves it must be remarked that the transition 
phenomenon is probably too unstable to be 
represented by a simple equation. 


11. PRAcTICAL APPLICATIONS 


Predictions of the maximum speed of airplanes 
based on wind tunnel tests involve a knowledge 
of the ‘‘scale effect’’ for the drag. Induced drag 


d turbulent friction (smooth plates). 


has practically no scale effect, also the drag of 
protruding parts follows fairly well the so-called 
quadratic law, at least beyond their critical R.N., 
which is connected with a more or less sudden 
drop of the drag coefficients. Thus the greatest 
uncertainty is that due to the scale effect on 
skin friction, i.e., the gradual drop of the friction 
coefficient with increasing R.N. The range of 
extrapolation is comparatively large in the case 
of modern fast airplanes. The comparison can 
be based on the R.N. of the wing defined in terms 
of the speed, the mean chord of the wing and the 
kinematic viscosity. Present atmospheric wind 
tunnels, with the exception of the Langley Field 
full scale tunnel, reach values of about R=2 
10°, the variable density tunnel about 3x 10°. 
The full scale values for large modern transport 
airplanes may reach 25 10°, somewhat smaller 
values corresponding to fast military planes. 
Hence there is a large change in the skin friction 
coefficient which influences the predicted speed 
very noticeably, especially because, in the case 
of modern clean airplanes, the skin friction may 
represent the major portion of the total drag. 
An exact calculation is made difficult by two 


facts: 
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(a) The separation of the skin friction and the 
drag of protruding parts is hardly possible. 

(b) The skin friction acts partly on the wings, 
partly on the fuselage and tail; these parts have 
different R.N. and so somewhat different scale 
effects. 

At the California Institute of Technology the 
following procedure, suggested by Clark B. 
Millikan, is in use: Two different values for the 
full scale drag coefficient are calculated, a first 
value based on the assumption that the total 
parasite drag follows the skin friction law for 
smooth plates, and a second value corresponding 
to the assumption that the profile drag of the 
wings follows the law mentioned, while for the 
remainder of the drag the scale effect can be 
neglected. It is believed that these values consti- 
tute limits between which the value of the actual 
drag should lie. 

As an example the following data are quoted, 
taken from the results of actual! wind tunnel and 
flight tests on an observation type, wire braced, 
gull wing monoplane. 


Model results scaled up without R.N. correction: 
Cp, = 0.0396, calculated Vinax = 180 m.p.h. 

Wing profile drag only corrected to full scale: 
Cp, =9.0364, calculated Vinax = 186 m.p.h. 

Total parasite drag considered as turbulent skin 
friction and corrected to full scale: 


Cp, =0.0268, calculated Vinax = 205 m.p.h. 


Flight tests on the actual airplane gave 
Vinax =196 m.p.h., i.e., almost exactly halfway 
between the two extrapolations. The relative 
distance from the two limits should, in general, 
be capable of estimation from the general 
“cleanness” of the airplane investigated. More 
refined methods of making such estimates have 
recently been developed by Millikan and will be 
presented in a paper to be published by him in 
the near future. 

Another important application of the skin 
friction theories is an approximate prediction of 
the boundary layer thickness or of the velocity 
distribution near the surface of plates or stream- 
line bodies. To be sure in the case of streamline 
bodies only a rough approximation can be 
obtained, because the ‘‘development of the 
boundary layer”’ is influenced by pressure distri- 


bution and curvature along the surface. How- 
ever, in many cases, an approximate calculation 
is of great value, for instance for an estimate of 
the drag of protruding parts in the boundary 
layer region, for an estimate of the influence of 
small roughness, or of the influence of the 
velocity field near the hull on the behavior of 
fins, rudders and propellers of airships. 

A simple formula for the boundary layer 
thickness can be obtained by comparing the two 
expressions for the coefficient of the local skin 
friction of a plate at the distance x from the 
leading edge, according to the Eqs. (33) and (38). 
We have the following two relations: 


logio (R.C;) (41) 
according to (38) and replacing in (33) R by R; 
log (R3(C,)'). (42) 


By subtracting the two expressions it is seen 
that the value R;/R,(C;)' is constant and 


R;=0.38R,(C;)! or 6=0.38x(C))!. (43) 


The value of C; as function of R, is tabulated 
on page 12; hence the value of 6 can easily be 
computed. The velocity distribution in the 
boundary layer is approximately given by the 
formula 


u = U[1—4.15(C,)! logis (6/y)]. (44) 


Let us calculate, as an example, the boundary 
layer thickness at the trailing edge of an airfoil; 
with chord of 3 m=9.8 ft. for the speed U=100 
m/sec. = 223.7 m.p.h. The corresponding R.N. 
(for ‘‘standard” air) = 20,000,000 =2 10"; C; 
= 0.00222, accordingly 6=5.4 cm. 


12. SKIN FRICTION OF BODIES OF REVOLUTION 
(AIRSHIPS) 


The application of the boundary layer theory 
to bodies of revolution is interesting for analyses 
and predictions of airship drags. The resistance 
of a bare airship hull is almost entirely due to 
skin friction, the pressures at the bow and the 
stern being practically balanced. Early wind 
tunnel tests, carried out at comparatively low 
R.N., appeared to be inconsistent with one 
another because the transition between laminar 
and turbulent flow was not taken into consider- 
ation. The airstreams in different wind tunnels 
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have very different degrees of turbulence and so 
the transition occurs in quite different ways in 
different wind tunnels. 

B. M. Jones and H. L. Dryden independently 
showed that such differences in airstream turbu- 
lence could account qualitatively for the wind 
tunnel discrepancies. In their considerations a 
figure of revolution was replaced by an ‘‘equiva- 
lent flat plate’ and the existing information on 
the flow along flat plates was thus applied to 
their problem. Shortly afterwards C. B. Millikan 
applied the boundary layer integral relation 
(analogous to Eq. (2), but taking the variation 
of pressure along the surface into account) to 
the case of axial flow about a figure of revolu- 
tion and gave expressions for the variation of 
boundary layer thickness and skin friction along 
such a body, assuming a parabolic velocity 
profile for the laminar portion of the boundary 
layer and a ‘“1/7th power” profile for the 
turbulent portion. The calculations were carried 
out in detail for the so-called N.P.L. models and 
curves of drag coefficient vs. R.N. for various 
degrees of airstream turbulence were obtained. 
These gave a satisfactory quantitative agreement 
with many of the experimental curves, so that 
the explanation of the latter’s inconsistencies 
was in this way successfully accomplished. 

Miss H. M. Lyon has carried the experimental 
side of the problem still further in an elaborate 
investigation of the correlation between the drag 
coefficient vs. R.N. curves for two airship hull 
models and the degree of turbulence produced 
artificially in the airstream. 

Freeman repeated Millikan’s calculations for 
the hull of the Akron dirigible and obtained a 
remarkably close agreement with his beautiful 
measurements on a 1/40th scale wind tunnel 
model, both as to drag coefficient and distribu- 
tion of boundary layer thickness. For the R.N. 
of his measurements the 1/7th power law used 
in Millikan’s theory was shown to be a very 
satisfactory approximation to reality. However, 
for the extrapolation of drag values to full scale 
and for full scale boundary layer calculations 
the 1/7th power law will almost certainly be a 
much less satisfactory approximation, so that it 
would be very desirable to repeat the analysis 
by using our improved knowledge of the laws of 
skin friction. 
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13. THE INFLUENCE OF ROUGHNESS 


The problem of the influence of surface rough- 
ness on the skin friction is rather complex. 
Although many experiments have been con- 
ducted concerning the friction in rough-walled 
pipes, considerably fewer tests are available on 
skin friction of rough plates in a free stream. 
In pipes with different diameters but with 
identical or geometrically similar surface condi- 
tions the coefficient of skin friction depends on 
the R.N. and the “relative roughness.’’ The 
relative roughness is the ratio between a length 
which is characteristic of the roughness and a 
length characteristic of the cross section (for 
instance, the hydraulic radius). It is not possible 
to define a priori which dimension is to be 
considered characteristic of the roughness, except 
that it must be the same for identical surfaces 
and for walls with geometrically similar kinds 
of roughness proportional to the scale of the 
protruding parts (for instance, to their mean 
height). Denoting this parameter by k and the 
radius by r, the friction coefficient will be 


c;= F(R, k/r). (45) 


The general form of this function is not known, 
in fact it is not even sure that one universal 
function covers the case of all kinds of rough 
walls. It is more probable that it depends on the 
shape of the roughness. Fortunately, the experi- 
mental evidence shows that the conditions are 
relatively simple for high R.N. and compara- 
tively large relative roughness. It appears that 
with increasing R.N. the skin friction coefficient 
becomes independent of the viscosity and in turn 
independent of the R.N. In this respect the 
behavior of rough surfaces is fundamentally 
different from that of a smooth surface, since in 
the case of the smooth surface the friction 
coefficient decreases indefinitely with increasing 
R.N. 

For the range in which the friction coefficient 
is independent of the R.N., a simple relation 
between the friction coefficient and the relative 
roughness can be found, based on the velocity 
distribution laws given in Section 6. We re- 
member that the ratio (velocity defect /friction 
velocity) was found to be a universal function 
of the ratio (distance/radius), and it has been 
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noted that this theorem holds also in the case 
of rough surfaces. Accordingly, denoting again 
the maximum velocity by U, the velocity at the 
distance y from the wall by uw, and the wall 
friction by 79, we write 


(46) 


The development for the velocity distribution 
starting from the wall will be different in the 
cases of smooth and rough walls. For smooth 
walls we assume that wu is a function of y, To, p 
and u. In the general case of a rough surface u 
will be a function of these quantities and also 
of the roughness parameter k. However, the 
fact that for large R.N. the skin friction coeffi- 
cient is independent of the viscosity indicates 
that the velocity distribution is also independent 
of v and hence u depends only on y, 79, p and k. 
Dimensional analysis shows that the only form 
possible for this functional relation is 


u=(ro/p)*g(y/k). (47) 


On the other hand, we obtained from the 
similarity of the turbulent flow pattern the 
relation 


u = log y) (48) 


so that due to (43) the development for the 
velocity distribution starting from the wall will be 


u = (to/p)*(const.+(1/x) log (y/k)). (49) 


Combining Eqs. (46) and (49) we obtain, by a 
reassumption similar to that applied in Section 8 


log (r/k)) (50) 
and 


log (r/R). (51) 


It is easy to show (compare Section 8), that the 
same equation holds for the coefficient referred 
to the mean velocity. 

The Eq. (51) was found by the author in 1930, 
by a somewhat different reasoning from that 
given above. The equation can be checked 
independently of the definition of k by tests in 
pipes or channels with walls of identical surface 
conditions. 

Let us write Eq. (51) in the form 


log k+(1/x) logr. (52) 


Then obviously, if we plot 1/c;} against log 7, we 


get parallel straight lines, each line corresponding 
to a given roughness. The slope of these lines is 
determined by the universal constant «x, which 
was encountered also in the case of smooth pipes 
and smooth plates. 

In Fig. 14 are plotted results of tests made by 
Fromm and Nikuradse. With one exception they 
agree with the foregoing conclusions. The ex- 
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Fic. 14. Plotting of 1/C;} for channels of different widths. 


ceptional case refers to a test in which a surface 
with extreme roughness (sawteeth of 1.5 mm 
height) was used. 

The vertical distance between any two such 
straight lines is determined by the difference of 
the terms (1/x) log k, i.e., by the ratio of the 
roughnesses. Thus, by choosing one definite 
roughness as a basis for comparison, it is possible 
to assign a definite value to any surface. Niku- 
radse produced geometrically similar rough sur- 
faces for all of his tests, making use of sand 
grains of various diameters. He found that by 
determining the roughness by Eq. (52), the 
roughnesses of the surfaces used by him are 
simply proportional to the diameter of the grains 
used. Thus we can tentatively use as a unit of 
roughness the grain diameter of an equivalent 
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Fic. 15. Plotting of 1/C;} against relative roughness. 


granular surface. Fig. 15 shows the results of the 
Fromm and Nikuradse tests combined into one 
straight line, which has the equation 


1/c;$=3.0+4.15 log (r/k). 


The next step is the application of Eq. (51) to 
the problem of skin friction of rough plates, in 
the same manner as we applied Eq. (33), deduced 
for the skin friction in smooth tubes, to the case 
of smooth plates. We must replace the radius 7 
in Eq. (51) by the boundary layer thickness. 
First, we notice that the ‘‘relative roughness” 
k/é which determines the friction is variable for 
constant k, it is large for small values of 6 (i.e., 
near the leading edge), and decreases with 
increasing distance from the leading edge. Hence 
the validity of such a calculation is restricted to 
cases in which the roughness & is so large that 
the assumption of a friction independent of the 
R.N. holds. 

We obtain, under these conditions, the equa- 
tion for the local friction coefficient: 


1/c3=5.8+4.15 log ((x/k)c;'). 


As might be expected, the friction depends on 
the ratio k/x (roughness relative to the length). 
In the following table are calculated a few 
corresponding values of c; and x/k: 


c,=0.0040 0.0035 0.0030 0.0025 
x/k=4xXi0® 8X10* 18X10* 53X10*% 10% 


(53) 


(54) 
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Fic. 16. Local friction coefficients for rough plates. 


These values give a fair estimate of the influence 
of roughness on the coefficient of skin friction 
for plates. Fig. 16 shows these values of c, for 
various roughnesses, compared with the curve 
for smooth plates. For a more exact calculation, 
it must be taken into account that in general 
the local skin friction will be a function both of 
the boundary layer R.N. and of the relative 
roughness. Experimental evidence for the range 
in which both the R.N. and the roughness have 
an influence on the friction is somewhat contra- 
dictory, even in the case of pipes. However, the 
interpretation of the experimental results is 
relatively easy in the case of the artificial 
granular roughness used by Nikuradse. He found: 

(a) that the influence of roughness is, for all 
practical purposes, null below a certain R.N., 
and that the limiting value of the R.N. increases 
with decreasing k. 

(b) that beyond this limit the friction coeffi- 
cient, for a given roughness, as a function of the 
R.N., first decreases, then slowly increases and 
becomes constant at a certain value of the R.N., 
which value also depends on the roughness. 

Considering the two limits (the limit at which 
the influence of the roughness appears, and the 
limit beyond which the friction coefficient is 
constant) Nikuradse’s tests show that both 
correspond to definite values of (70/p)!k/v. This 
can easily be understood: Let us remember that 
in the case of a smooth wall, we have in the 
immediate neighborhood of the wall a laminar 
region, whose outer limit or “‘the thickness of the 
laminar layer” 6; is proportional to v/(to/p)}; 
thus the parameter (70/p)'k/v is proportional to 
k/6, or to the ratio between the roughness and 
the thickness of the laminar sublayer. 
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According to the tests, the influence of the 
roughness starts with (70/p)'k/v=3, and the 
friction coefficient becomes constant for (70/p)*k/ 
vy >60. We remember that 6;211.5y/(70/p)! gives 
about the ‘‘thickness” of the laminar sublayer. 
Thus it appears that the influence of roughness 
becomes noticeable if the height of the roughness 
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TABLE III. 
x k x k 
0.075 m 0.0230 mm 2m 0.0306 mm 
0.1 0.0237 2.5 0.0311 
0.25 0.0257 3 0.0316 
0.5 0.0276 10 0.0344 
1.0 0.0290 50 0.0385 
> 0.0299 100 0.0400 
200 0.0418 


has the same order of magnitude as the laminar 
sublayer (k/6;>4) and that the viscosity has no 
more influence if the roughness is large in 
comparison with this thickness (x/5,>6). 

Accepting this evidence for the influence of 
roughness with the limiting value (70/p)!k/v=3, 
we can carry out a quite interesting calculation 
indicating the magnitude of ‘‘allowable rough- 
ness’’ in practical cases—‘‘allowable”’ in the 
sense that roughness below this limit does not 
influence the friction of a plate. 

According to this condition the limiting value 
of k is given by 


k (55) 


Now, for smooth plates c; is a function of the 
R, and is calculated on page 12. Using these 
tabulated values, we obtain for the case of an 
airfoil with 3 m chord and 100 m/sec. speed 
(x=distance from the leading edge), Table II. 


TABLE II. 

x k x k 
0.075 m 0.0100 mm 1.5 m 0.0129 mm 
0.1 0.0103 2.0 0.0131 
0.25 0.0111 2.8 0.0133 
0.5 0.0118 3.0 0.0135 
1.0 0.0125 


Similarly, for an airship of 200 m length and 
speed of 40 m/sec., we have Table III. 

It is interesting to notice that recent experi- 
ments of the N.A.C.A. give the same order of 
magnitude for the allowable roughness of an 
airfoil. 

It is evident that a certain small roughness is 
more harmful near the leading edge and less 


harmful at the rear portion of the surface. This 
has been proved experimentally by G. Kempf 
and corresponds to common experience. 

Another consequence which can be drawn from 
Eq. (55) is the fact that the allowable rough- 
ness k depends chiefly on the speed and only 
to a smaller degree on the R.N. This enhances 
the reliability of wind tunnel tests, since it shows 
that in model scale-reduction it is not necessary 
to reduce k in the ratio of the length scale, but 
only in the ratio in which c;} is increased. 

Unfortunately the evidence of tests on pipes 
with different roughnesses produced by various 
paints or by different manufacturing methods is 
not so consistent as are Nikuradse’s tests on 
artificially roughened pipes. Some of these other 
tests are very exact, hence it is probable that 
Nikuradse’s conclusion as to the existence of a 
universal ‘‘roughness function’? does not cover 
all kinds of roughness. It is probably that the 
one parameter & is not sufficient in itself to 
characterize the surface; at least one more 
parameter, such as the mean distance between 
protruding parts, must enter into the picture. 
At the present time it is not possible to compute 
the skin friction for a plate with a given rough- 
ness. 

The author believes that the main problems 
of future research in the field of skin friction are: 

(a) to work out a method of determining a 
reliable “roughness scale” for different kinds of 
materials used in practice, 

(b) to investigate, by systematic experiments 
on rough surfaces, the transition range, i.e., 
the range in which both the R.N. and the 
relative roughness influence the friction. 
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Air Force Measurements on Bodies Moving Through Still Air 
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HE six aerodynamic reactions, consisting, 

in the most general case, of three forces 
and three moments which are sustained by a 
body moving through a fluid or by a body held 
in a fluid stream, may be classified on the basis 
of their origin into two general groups, those 
which are due essentially to surface pressure and 
those arising principally from surface friction. 
Wing lift and pitching moment are common re- 
actions of the first group, while profile drag and 
yawing moments are examples of the second. 
Generally speaking, the pressural reactions on a 
body of easy form and low incidence are very 
nearly the reactions associated with potential 
flow and hence are quite closely proportional to 
the square of the velocity. They are, therefore, 
essentially free of scale effects and possess a 
certain degree of stability against the various 
disturbances which are found in the air streams 
of wind tunnels. In such cases, the pressural 
reactions which are measured on a model in a 
wind tunnel, when multiplied by the square of 
the ratios of scale and speed, are the reactions 
on the full scale body. There are other cases, 
however, where the pressural reactions as meas- 
ured in the wind tunnel have no simple relation- 
ship to the full scale values. 

When the burble attitude! or the critical 
bluntness? of a body at which the flow breaks 
away from the surface is approached, the 
pressural reactions become unstable and the flow 
regime becomes sensitive to various disturbances 
such as the turbulence initially present in a 
tunnel stream.’ In these critical cases the close 


1A wind burbles, for example, at —6° or —8° angle of 
attack referred to the null line and again at +14° or 
+16°. 

2 The critical bluntness of a spindle form or airship hull 
is somewhat uncertain but appears to be near a fineness 
ratio of two. 

3 For the effect of tunnel turbulence on wing burbling, 
for example, see Effect of Turbulence, by Millikan and 
Kline, Aircraft Engineering, August (1933). 
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measurement of the pressural reactions on a 
body in a wind tunnel, regardless of its size, 
are either of doubtful value or are distinctly de- 
ceptive and cannot be safely used as a basis 
for the study of important dependent phenom- 
ena such as spinning or lateral control at the 
stall. 

The frictional reactions present a more serious 
problem because of their sensitivity to turbulence 
and other tunnel disturbances not only at the 
burble attitudes and critical fineness ratios, as in 
the case of the pressural reactions, but at all 
other attitudes and fineness ratios as well. The 
general dependence of the frictional reactions at 
a given Reyrfolds number upon tunnel dis- 
turbances is due essentially to the effect of such 
disturbances upon the position of the transition 
region which separates the laminar and turbulent 
phases in the boundary layer. The most im- 
portant imperfections of the wind-tunnel air 
stream in this connection are the presence of 
pressure gradients alongstream and initial stream 
turbulence. The effect of the first is to retard the 
transition and that of the second to advance it, 
the latter effect normally being much the more 
important. Other disturbances, all of which tend 
to advance the transition, although to a lesser 
extent and hence to increase the drag coefficient 
at a given Reynolds number, are shrill or ex- 
plosive sounds, fluttering of the body or of the 
air stream pressure and high frequency vibra- 
tions of the body surface. Fortunately, these 
latter disturbances can all be practically elimi- 
nated without difficulty from wind-tunnel air 
streams. 

The uncertainty in the position of the transi- 
tion from the laminar to the turbulent boundary 
layer phase, as it is effected by these disturbing 


factors, is quite different for different bodies and 


at different Reynolds numbers. The transition 
position for flat plates, expressed in terms of 
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the transition Reynolds number,‘ may occur, 
for example, at any value between 10° and 
510° depending upon the kind, the amount 
and the number of the disturbing factors.® 
A plate nine inches wide in an air stream whose 
velocity is 100 feet per second may have a 
boundary layer, therefore, which is all laminar 
or 20 percent laminar and 80 percent turbulent, 
depending upon the wind tunnel and the model 
mounting, the relative drag coefficients being for 
the two cases about one to two. 

In the case of spindle forms, such as airships, 
which have a surface point of minimum pressure 
near the bow, followed by a region of rising 
pressures and usually by a second point of 
minimum pressure near the stern, the position 
of the transition region is not only very sensitive 
to disturbing factors, over the common wind- 
tunnel range of Reynolds numbers but also tends 
to be unstable. The instability at these Reynolds 
numbers is due to a tendency of the transition 
region to leap upstream across the midship area 
of rising pressures to the forward minimum 
pressure point when turbulence appears or in- 
creases in the air stream. Such a shift in the 
transition region may double the resistance 
coefficient of the body at a given Reynolds 
number. Under such circumstances the accurate 
forecasting of full scale resistance coefficients 
for airships from measurements made in wind 
tunnels is quite hopeless. 

For these reasons the initial turbulence in 
tunnel air streams has come to be recognized by 
present day experimental aerodynamics as an 
intractable if not an insolvable defect in the 
wind-tunnel method. Notwithstanding notable 
advances which have been made, turbulence 
still cannot be successfully measured or pre- 
vented, neither can it be completely controlled 
or defined. Its existence in varying degrees in 
different tunnels continues to be the chief cause 


4 The transition Reynolds number is defined as VL; /» 
where L;, is the distance from the bow, or the leading 
edge, to the transition region between the two boundary 
layer phases. 

5 Some important specific evaluations of the transition 
Reynolds number are 1.6 10° by Elias, N.A.C.A. Tech. 
Memo. No. 614; 3.210° by Hansen, N.A.C.A. Tech. 
Memo. No. 585; and 4.8X10° by Burgers and van der 
Hegge Zignen, Delft (1924). 
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of a variety of inter-tunnel disagreements and 
inconsistencies. 

In order to evade the whole turbulence prob- 
lem, the Massachusetts Institute of Technology 
became interested early in 1930 in conducting 
tests in calm, non-turbulent air. Under the 
sponsorship of the late President Stratton, who 
was chiefly interested in the project, plans were 
laid for the use of the proposed model basin, 
the completion of which was then in prospect. 
It was considered advisable, however, to conduct 
some preliminary studies in the interim. After 
several methods were discussed the following 
apparatus was agreed upon for these preliminary 
investigations. 

There exists at the Institute a large room, or 
rather open well, in the Steam Engineering 
Laboratory, the size and shape of which may be 
seen in Fig. 1. Although the air here is more or 
less disturbed during the day due to convection 
currents from hot steam engines on the lower 
levels, drafts from open windows, etc., the air 
is very quiet and the temperature quite steady 
late at night after the building has been shut up 
for some hours. A number of observations of 
smoke from titanium tetrachloride have shown 
that from about 11 P.M. until 6 A.M. the maximum 
velocity of the air in any direction is about 0.2 
of a foot per second and that it is often so still 
that no movement can be detected. This motion 
is regarded as negligible in comparison with the 
speed at which the tests are run. An attempt to 
secure a quantitative measurement of the tur- 
bulence existing in this body of air with the 
methods at present available did not seem likely 
to produce results of much value. Although it is 
probable that some very slight degree of large 
scale turbulence is present,® the air appears to 
be as near to the condition of complete rest as it 
is practical to come with so large a mass of air. 

The method agreed upon for the preliminary 
tests was to move the models through the still 
air of this large room in a straight line with as 
nearly a uniform velocity as possible and to 
measure the speed, the model lift and the 
drag simultaneously. The models were to be 


6 Tetrachloride smoke, for example, moves in a very 
leisurely manner and seems to have very little tendency to 
form eddies or whorls. 
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Fic. 1. The site of the coasting apparatus with principal dimensions. 


mounted on a suitable carriage running on two 
long and very taut inclined wires; the model 
and carriage were to be launched at nearly the 
terminal velocity, allowed to coast down the 
wires and finally brought to rest by a suitable 
retarding gear. The absolute speed was to be 
measured with an electric tape chronograph, 
while the lift and drag forces were to be recorded 
as continuous records on a smoked glass chart 
borne and rotated by the carriage, and the two 
records correlated. The drag mechanism, at 
least, was to be in dynamic balance and quite 
independent of the lift mechanism. Both were 
to be oriented accurately to each other and to the 
wire tracks. 

The final apparatus which was built to these 
general specifications is shown in Fig. 1 which 
outlines the wire site, and in Fig. 2 which gives 
a close view of the travelling balance. The main 
span is of two 1/8-inch high strength alloy steel 
wires, one placed 8 inches above the other and in 
horizontal projection is 213 feet long. The wires 
are wedged into steel blocks at the ends in such 
a way that the fastening will develop the full 
strength of the wire. At the lower end of the span 
the wire anchor blocks are fastened to an elevator 
which is supported on a structural steel tower. 


This elevator is suitably guided on the tower 
and is provided with a motor, counterweight, 
etc., so that it may be driven up or down the 
tower under its own power by an operator who 
rides with it. By moving the elevator up or 
down the slope of the span wires may be changed 
from approximately 1° minimum to 10° maxi- 
mum. Wire tautness, sufficient to reduce the 
bare wire sag to about one inch, is maintained by 
two 850-lb. double reeved weights, which are 
attached to the upper end of the span wires by 
means of a suitable sheave arrangement. 


Fic. 2. The travelling balance ready to run. 
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Fic. 3. The drag system of the travelling balance. 


The travelling balance, Fig. 2 and shown 
diagrammatically in Figs. 3 and 4, performs the 
following functions: It carries the model; it 
resolves the aerodynamic reaction on the model 
accurately into its lift and drag components and 
records their values; it sends out signals at 
exactly equal distance intervals along the span 
and simultaneously makes a tally mark on the 
smoked glass disk by which the records of speed 
and air forces can be correlated. The drag is 
measured by the displacement of the drag truss 
A, Fig. 3, which is supported on four rolling 
sectors B against the drag spring C and is 
recorded on the smoked chart D through a fine 
wire and a pointer. A vital part of the drag 
mechanism is the dynamic balance, marked E in 
the figure, whose function is to suppress all 
effects of accelerations in the direction of flight. 
It consists essentially of a mass supported on 
two rolling sectors and closely guided so that its 
path and that of the drag truss are accurately 
parallel. This mass carries a roller which works 
in one end of a forked lever, F; the other end of 
the forked lever works on a second roller which 
is fastened to the drag truss. By this arrange- 
ment the acceleration force of the total mass of 
the drag truss system acting at the truss roller 
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Fic. 4. The lift system of the travelling balance. 


and the corresponding force of the mass // 
acting at roller L can be made to give equal and 
opposite moments about the point G by a simple 
adjustment of the arm length GL. When roller L 
is so adjusted by means of the micrometer screw 
provided for the purpose, the drag system is 
free of all effects of translational accelerations in 
the direction of the drag. No attempt was made 
to balance the moments due to angular accelera- 
tions, since such accelerations are negligibly 
small. 

The lift is measured by the displacement of the 
lift system, Fig. 4, against the lift spring. 
Member K is carried by the drag truss but is free 
to move perpendicular to it, the displacement 
being recorded by a pointer system on the same 
smoked chart that records the drag truss dis- 
placement. 

Particular attention has been given to the 
reduction of the pitching of the travelling 
balance to a value that will not impair the 
precision of the drag measurements through an 
error in the force resolution. For this reason the 
carriage has been made quite long and is de- 
signed to carry its models below the center 
section so that the resultant forces on lifting 
models will pass near the center of gravity of 
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the carriage and induce a minimum of pitching. 

The entire instrument, Fig. 2, is constructed 
from rolled duralumin shapes and high strength 
aluminum alloy castings. It is very light, weigh- 
ing only 14.5 pounds with housing, yet is 
sufficiently rugged and stiff to maintain its 
internal alignments and adjustments in normal 
usage and, in the case of an emergency, to with- 
stand without serious injury a force of 200 or 
250 pounds, corresponding to a deceleration of 
15 or 18g, applied at the bow. All working parts 
are fully shielded from the air by a transparent 
celluloid housing streamlined and proportioned 
to render the interference effects on the model 
negligibly small. 

The practical speed range of the travelling 
balance is from 5 to 40 feet per second and is 
controlled by adjusting both the launching speed 
and the wire slope and, when necessary, by 
adding aerodynamic drag to the carriage. No 
difficulty is found in adjusting the speed so 
that it is uniform within +0.1 foot per second 
over the central span region. The small wire sag 
gives a slight decline in the gravity towing com- 
ponent along the span which introduces a small 
but very useful variation in the velocity, since 
it makes it possible by launching slightly under 
the terminal speed to obtain a speed maximum, 
and hence speeds of nearly constant magnitude, 
in the mid-span region. 

The method of measuring speeds is simple and 
accurate. As the travelling balance advances, a 
trolley attached to the balance frame and held 
in contact with the upper wire by spring pressure, 
rotates the glass chart, through a suitable gear 
train, at the rate of one revolution per run. It 
also rotates a Bakelite cam with a metal insert 
at the rate of twenty revolutions per run. When 
the cam follower makes contact with the metal 
insert, the chronograph circuit is closed and a 
mark is made on the chronograph tape. Simul- 
taneously, a scriber on the cam follower shaft 
makes a sharp break in its record line on the 
smoked chart, Fig. 5. Since the tape speed is 
constant—three, six or twelve inches per second 
depending upon the need—the length of tape 
between any two adjacent marks is a measure of 
the average velocity for the corresponding in- 
terval of span. The part of the smoked glass 
record which was made in passing the particular 


bo 


Fic. 5. A typical smoked glass record. 


interval can be easily identified. Under normal 
conditions when the speed is nearly constant, 
the speed precision is approximately +0.5 per- 
cent. The chronograph circuit’? is shown in 
Fig. 6. 

The gear for launching the travelling balance 
at various speeds is shown diagrammatically in 
Fig. 7. The energy for launching is supplied by 
a falling 500-lb. weight which transmits its 
energy to the balance through a system of 
pulleys and cables, as shown in the figure. The 
launching acceleration is about 4g which yields 
a maximum speed of about 40 feet per second 
when the catapult travels its maximum distance. 
For intermediate launching speeds the catapult 
is pulled back to an intermediate position meas- 
ured on a suitable fixed scale. 

The arresting mechanism, shown diagram- 
matically in Fig. 8, is divided into two parts, 

7 The chronograph circuit arrangement is as follows: The 
two insulated span wires are connected into the grid 
circuit of a vacuum tube amplifier so that the ‘‘make”’ 
puts a heavy negative bias on the grid, thus allowing 
current to flow in the plate circuit and operate the chrono- 
graph pen. The arm which carries this pen is provided with 
an interrupter which breaks the circuit as soon as the pen 
moves. The plate circuit is then dead and, since the 
travelling balance no longer holds the grid circuit closed 
after passing over the station, the grid takes a positive 
potential from battery B, Fig. 6, and further operation of 
the plate circuit is prevented until the next station is 
reached. In this way the pen is prevented from making 
more than one jab for each station. 
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the service gear and the emergency gear. The 
Vert rere service gear consists of two lengths of sash cord 
~ SSOCEGREREE which cross the flight path and, after running 
calcite oe through a series of pulleys finally wind on a 
hg ee small drum which is provided with a band 

Nt brake. The band is anchored at one end and 
) THYRATON F627 pulled tight by a hanging weight at the other. 
. CHRONOGRAPH The anchorage is adjustable so that the angle of 
| wrap may be varied by the operator who can 
ly thus control the “hardness” of the gear. The 
2 emergency gear is behind the service gear and 

| consists simply of three airplane shock cords, as 
shown. 

A typical smoked glass record, with certain 
information added to indicate the method by 
which the forces are evaluated, is shown in 
Fig. 5. The records are regularly read with the 
aid of a low power microscope mounted on a 
swinging arm and so located with respect to the 
record that the path described by the cross hairs 
coincides with the arcs made by the original lift 
or drag scribing point. A suitable scale is pro- 
vided by which the position of a record line 
which runs between two calibration circles can 
be quickly and accurately read. 

The precision of the drag readings under the 
conditions of static calibration is +1 gram for 
drags under 50 grams, and +2 grams for drags 
between 50 and 500 grams. The lift readings are 
reliable under the same conditions to +2 grams 
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Fic. 6. The chronograph circuit. 


CATAPULT LOADING 


for loads under 100 grams, +3 grams for loads 

between 100 and 500 grams and +5 grams for 
loads between 500 and 1500 grams. Calibration ' 
test runs were made to determine the minimum 


distance between model and carriage necessary 
to avoid measurable interference. 

As a part of the calibration of the apparatus | 
investigations were carried through on two sets 
of models, a series of circular disks of varying 
diameters and a series of equiarea rectangles of 
varying aspect ratios. All were held normal to 
the wind. Since the drag of such bodies is entirely 
pressural in origin and since the point of break- 
away is fixed at the body edges with great 
stability, the drag must be nearly independent 
of initial turbulence in the air. The wind tunnel 
and the calm air drags of such bodies should, 
therefore, agree. It will be seen that the experi- 
mental results show an agreement within fair 


Fic. 8. Diagram of arresting gear. precision. 
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The series of circular disks, consisting of five 
disks ranging in diameter from 2 to 13 inches, 
was made from 1/16” thick 17st aluminum 
alloy machined circular to the chosen diameter 
to a tolerance of +0.005’’. The edges were made 
square with the thickness of the outer portion 
reduced so that all disks were in effect geo- 
metrically similar with ¢/d=208. For con- 
venience, the diameters were chosen so that the 
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intervals between the logarithms of the diameters 
would be equal. This series was tested in the 
five-foot M. I. T. wind tunnel and on the 
coaster. The results are shown in Fig. 9. Agree- 
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Fic. 9. Drag of circular disks in wind tunnel and in 
still air. 


ment is excellent over a considerable range, well 
within +3 percent for Reynolds numbers be- 
tween 80,000 and 225,000, although at very low 
Reynolds numbers, the precision of both wind 
tunnel and coaster begins to decline, due to 
the smallness of the forces. In Fig. 10 the 
results of other studies are also given for com- 
parison. 

The second series consisted of nine rectangles 
all having the area 51.8 square inches but varying 
in aspect ratio from 1 to 40. The calm air drag 
coefficients are plotted in Fig. 11 along with 
the coefficients which were measured in a tunnel 
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Fic. 10. Drag of circular disks. Comparison of still air 
results with values reported by others. 


by Prandtl® and those found from drop tests 
by Eiffel. The agreement is fair, varying from 
complete agreement at the aspect ratio 4 to a 
departure of about 5 percent at the aspect ratio 
15 and of nearly 10 percent at the ratio 1, the 
wind-tunnel values being low. It is noteworthy 
that the drag coefficients appear to approach the 
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Fic. 11. Drag of rectangular plates. Comparison of coaster 
results with those reported by others. 


value 2.0 asymptotically with the increase of 
aspect ratios toward infinite values, and that 
the gradient of the coefficient with aspect ratio 
at the ratio unity is now found to be zero as it 
should be. 


5L. Prandtl, Ergebnisse der Aerodynamischen Versuch- 


sanstalt zu Gottingen. 
Eiffel, The Resistance of the Air and Aviation. 


Translated by J. C. Hunsaker. 
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Aerodynamics of the Wing Slot 


Max M. Munk, Catholic University of America 
(Received August 15, 1933) 


HE large maximum lift coefficient of a wing 

equipped with a slot parallel to its leading 
edge is a discovery of sufficient importance to 
warrant more study and research than it has 
hitherto received. The addition of the slot gives 
the wing a larger lift than one would expect, and 
further analysis may lead to more efficient slot 
arrangements. A rational explanation of what 
happens will not only satisfy the intellectual 
curiosity but will invite suggestions for further 
research which, if properly conducted, is bound 
to lead to valuable improvements in airplane 
design. 

Aspect ratio which is so important in con- 
nection with the drag and the lift at a specific 
angle of attack has little, if any, relation to the 
large lift in the case of slotted wings. This fact 
simplifies the discussion to that of flow about a 
wing in two dimensions. For all ordinary wing 
sections the lift coefficient may be considered as 
increasing linearly with the effective angle of 
attack by about 0.1 per degree until the angle 
of attack of the chord has reached about 10°. 
With a further increase in angle of attack, there 
results a rapid change in the flow pattern, 
accompanied by an increase in the drag, the lift 
remaining the same or even decreasing. 

It is of importance to see how this brief ac- 
count of the facts underlying the theory of the 
unslotted wing must be altered to apply in the 
case of slotted wings. If the unslotted wing 
theory is applied to each component wing in- 
dividually, first to the narrow front wing, here- 
after called the slat, and then to the main wing 
in rear of the slat, some progress will be made 
towards explaining the larger lift caused by the 
slat. It will be found, however, that the total 
effect computed by this method is not large 
enough to agree with experiment. 

First, let us consider whether it is possible for 
a slotted wing, when considered as a whole, to 
have a maximum lift coefficient larger than the 
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maximum lift coefficient of its components. In 
other words, if both the slat and the main wing 
have maximum lift coefficients of say 1.5, based 
on their separate areas, the question is: can the 
composite wing have a maximum lift coefficient 
based on the total area of say 2.5? Obviously it 
cannot as long as all three lift coefficients, namely 
that of the slat, that of the main wing and that 
of the composite wing, are computed by the use 
of the same velocity, namely the velocity of 
flight. This raises the question whether use of 
another velocity, in the case of the slat, would 
not comply better with the wing-section theory, 
and whether the more correct velocities are suf- 
ficient to explain the larger lift. 

The use of the local velocity, rather than the 
velocity of flight, is common practice in the 
older papers on biplane theory. The air velocity 
at each element of a composite wing should not 
be considered, according to the biplane theory, 
to be the velocity of flight but rather that 
which includes the interference components due 
to the adjacent wing or wings. This method 
should also be valid for the components of a com- 
posite wing. 

First, consider the air force of the main wing 
which experiences the interference of the slat. 
It can be seen at once that the interference 
effect on the main wing cannot be large, since 
the narrow slat is unable to set up velocities 
extending far enough to influence materially the 
much larger main wing, and since the local 
velocities in the rear of the slat or of any wing 
are low. Hence, in the light of the wing-section 
theory, we cannot expect an appreciable increase 
of either the lift or the formal lift coefficient 
(computed with the velocity of flight) of the 
main wing. This is so obvious that it is not 
necessary to support it by mathematical detail. 
The application of the wing-section theory to the 
main wing disappoints the investigator as a 
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rational explanation of the entire lift increase as 
found experimentally. 

The experimental results obtained by means of 
pressure distribution measurements by A. Or- 
merod,! for example, show likewise that ‘‘the lift 
coefficient of the main wing alone, operating be- 
hind the slat, never exceeds the maximum lift 
coefficient of the unslotted section’ and that 
the slat bears the burden of the increase in lift. 
This confirms our contention. 

The explanation of the slat lift, however, 
requires more study. In view of the relatively 
large chord of the main wing, the entire slat can 
be considered to be in the immediate neighbor- 
hood of the main wing, although the reverse does 
not hold. The slat is, moreover, located near the 
front or leading edge of the main wing where the 
local velocities which are set up by the main 
wing are known to be large. Let us form an 
estimate of the magnitude of these local velocities 
and then compare them with the interference 
velocities which are necessary to explain the 
large lift of the slat. 

Let the velocity in front of a wing with an 
angle of attack of 30° be estimated for a position 
about 5 percent of the chord in front of the 
leading edge. This would correspond to the 
proportions in Ormerod’s tests. Since an estimate 
of the velocity is sufficient for this explanation, 
one can proceed with a minimum of mathematical 
analysis. The amount of work involved is 
lessened not only by segregating the dominant 
physical effects from the non-dominant ones but 
also by segregating the dominant mathematical 
terms from the non-dominant terms. This 
approximate solution will enable one to compare 
the local and interference velocities and to 
explain the large lift on the slat; further sepa- 
ration would be wasted. 

The analysis proceeds by substituting the 
nearest ellipse for the wing section and by sub- 
stituting for the ellipse a straight line connecting 


' the two foci. The component of the flow velocity 


normal to the straight line should not be notice- 
ably different from the normal component of the 
velocity of flight, that is, V sin a. Now, in the 
plane of complex numbers, let the straight line 
extend between the two points +7. The velocity 


' A, Ormerod, R. and M. No. 1477, p. 4, line 3. 
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function of the flow set up by a lateral motion 
of this straight line is represented by the velocity 
function, 


F,=V sin a(1—2/(1+2%)!), 


where z denotes the complex variable repre- 
senting the points of the plane. This function 
gives zero velocity at infinity, and between the 
foci points on the straight line, where z is +7, 
the second term in the bracket becomes purely 
imaginary and hence the real component becomes 
constant as required. Since the velocity dis- 
tribution near the center, z=0, does not vary 
much, the same flow will also serve approxi- 
mately for elongated ellipses with the ends of 
the straight line as foci. The value of z repre- 
senting the region of interest is substituted in 
the expression for the velocity function. Now 
the distance of the center of the slat from the 
leading edge or from the focus is about 5 percent 
of the chord, since the distance between the 
leading edge and the adjacent focus can be 
neglected in the present case of a very elongated 
ellipse. We have, therefore, to substitute s=1.1 
into our expression for the velocity function and 
2 for the length of the straight line representing 
the chord. This gives an interference velocity of 
the order of magnitude. 


V sin a(1—1.1i/(1 —1.12)!) =0.76V sin a=0.38V, 


since the angle of attack is about 30°. 

This, however, is only a portion of the inter- 
ference, since we have not yet taken into account 
the circulation flow of the wing near the leading 
edge. The circulation interference, under ordi- 
nary conditions, is about equal to the inter- 
ference of the transverse flow. In the present 
case, the circulation could correspond, at best, 
to an angle of attack of say 10°, because the lift 
of the main wing corresponds to that angle of 
attack. The total interference would then be of 
the order of magnitude of half the velocity of 
flight. It would be larger in rear of the slat, 
because the trailing edge is closer to the main 
wing and smaller at the leading edge of the slat; 
perhaps the resultant effect would correspond 
to a value greater than the interference at the 
center of the slat. The difference, however, would 
be small. It should be noted that the interference 
just mentioned must be superposed on the com- 
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ponent of the velocity of flight parallel to the 
wing section. This component, which is slightly 
smaller than the velocity of flight, is not parallel 
to the total interference velocity. Taking these 
facts into consideration, we are forced to the 
conclusion that the local velocity of the slat will, 
in no case, exceed one and one-half times the 
velocity of flight. Furthermore, the slat increases 
the apparent lift coefficient of the wing by about 
0.7 when its chord is 14.5 percent of the chord of 
the main wing. Hence the lift coefficient referred 
to the area of the slat would be increased by 
0.7/0.145, or approximately by 5, or say from 
1.6 to 6.6, which is slightly more than a four-fold 
increase over its original value. Also, since the 
lift varies as the square of the velocity, the local 
velocity must be about twice the velocity of 
flight. The former considerations of local velocity 
explain, at best, half the lift increase. 

There must then be a second effect not associ- 
ated with single wings in free air. It cannot be a 
change of the direction or the curvature of the 
local air flow since these two influence only the 
lift at a given angle of attack and not the 
maximum lift obtainable. There remains then 
only the local rate of variation of the velocity, 
the gradient of its magnitude and the difference 
of direction of the velocity between adjacent 
points. As a rule, the rear edge of the slat is 
located in a region of rather high local or inter- 
ference velocity, while the leading edge of the 
slat is in a region of much smaller local velocity. 
The air flows through the slot in a convergent 
manner. It is as if the slat were located in a tube 
of flow of narrowing cross section. It is, however, 
well known that convergent air motion en- 
courages a smooth and orderly air flow. In front 
of an obstacle much less streamlining is required 
than in the rear, because the air is convergent 
in front and divergent in the rear. For the same 
reason, a venturi tube has to be longer aft the 
smallest cross section than in the front. It 
therefore stands to reason that the convergence 
of the local air flow along the slot is most 
probably responsible for the lift increase other- 
wise unexplained. The burbling of the fluid in 
the case of a wing with a slat is delayed in 
consequence of the convergence of the local flow. 
This is as far as reasoning will permit us to go. 
Whether or not the explanation is a correct one 
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can only be proved by experiments. Wing 
sections should be investigated in convergent 
flow in this connection to see whether their 
maximum lift coefficient can be increased. 
There is a second very interesting test sug- 
gested by the preceding analysis. A slotted wing 
is considered as the combination of a small wing 
and a large one, the small one being in the local 
field of the large one and under the influence of 
the increased velocity and convergence created 
by the large wing. This is similar to a multiple 
venturi tube where a small venturi is placed at 
the narrow cross section of the larger one in 
order to utilize its large local velocity. This 
method of increasing the vacuum has_ been 
carried further in this case, by placing a third 
tiny venturi into the second one at its smallest 
cross section. The analogy to the wing is obvious. 
It would be very instructive and interesting to 
place an even smaller slat near the front of the 
present slat in a local region of large velocity and 
convergence and to observe the resultant lift. 
As with the venturi tube, the principle cannot, 
of course, be carried too far. We cannot put one 
hundred venturi tubes one into the other, not 
only because that is not practical but also 
because the effect would, in the end, be lost. 
The pressure differences would assume magni- 
tudes large enough to obtain compressibility 
effects. The velocities would approach the 
velocity of sound where the ordinary laws of a 
perfect and incompressible fluid no longer hold. 
This analysis may be concluded with a few 
words about a quite different theory or explana- 
tion of the wing slot effect. L. Prandtl and A. 
Betz? explain the lift increase of the slot by its 
effect on the boundary layer. They say that the 
air stream flowing through the slot lifts or raises 
up the boundary layer. In other words, an 
attempt is made to explain the slot effect by di- 
rect air-friction influences, the argument being 
that by subdividing the wing the friction effect is 
divided. The main wing remains subject to the 
effects of its own surface friction but is no 
longer under the influence of the friction along 
the slat surface. On the contrary, the friction 
layer, which means the layer of air strongly 


2 L. Prandtl and A. Betz, Ergebnisse der Aerodynamtischen 
Versuchsanstalt zu Gottingen, Lieferung 3, p. 9, line 17. 
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affected by the surface friction of the slat, 
passes by the main wing at a safe distance 
without unfavorably affecting the air flow. 

This explanation seems unsatisfactory in 
several respects. According to it, an unslotted 
wing would not have less desirable charac- 
teristics than the same wing behind a slat, hence 
the latter wing should not have any larger lift 
than an unslotted wing, as indeed it has not, 
and we are exactly where we were before. Fur- 
thermore, since the effect of the slot on the main 
wing and not the effect on the slat is discussed, 
the explanation would imply that the main wing 
bears the burden of the lift increase, whereas the 
experiment shows that the slat does. Further, 
if the effect is a pure and direct friction effect, 
the same result for both cases would be obtained 
by diminishing the influence of the friction, for 
instance, by increasing the model size. Experi- 
ence shows that the maximum lift coefficient is 
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not increased by an increase of the Reynolds 
number but rather decreased in keeping with a 
much broader principle. 

Finally, it would seem that an explanation 
based on the application of a generalized wing 
theory and applied to the components of the 
composite wing is preferable to an explanation 
based on direct friction effects. We do not know 
enough about air friction. The fundamental 
problem of the motion of slightly viscous fluids 
has not yet been solved. Any explanation based 
entirely on air friction can only explain one 
unknown quantity by another unknown quan- 
tity. This is not only logically unsatisfactory but 
also leaves the alleged explanation useless for the 
further advancement of science. A reduction of 
the problem to an application of a wing theory 
applied to each individual wing will lead directly 
to a chain of constructive suggestions for further 
theoretical and experimental research. 
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NE of the greatest obstacles to progress in 
meteorology is undoubtedly to be found in 
the wide gulf between the mathematical theory 
on one hand and the applied science, weather- 
map analysis and forecasting, on the other. In 
aerodynamics, it is possible to advance a theory 
and then to proceed to a set of specially arranged 
experiments from which the validity of the theory 
may be conclusively established. In meteoro- 
logical research, experimental verification is 
usually impossible and one is therefore forced to 
study the full scale experiments provided by the 
atmosphere. In this study the chief difficulty 
lies in the number of variables involved in even 
the simplest phenomenon. 

To eliminate all but one of these variables, two 
methods are possible. Utilizing a great number of 
observations, one may group the observed effect 
with respect to one selected variable, assuming 
the random distribution of the other elements to 
be the same in each group. This might be true 
if assurance could be had that the variables 
really are mutually independent but this, un- 
fortunately, is seldom the case. 

The other method consists in a study of 
individual cases, a procedure which has gained 
increased favor with the advent and ascendance 
of the Norwegian Meteorological School. The 
chief drawback to this method consists in the 
obvious opportunities it offers for a biased 
selection of evidence in favor of a preconceived 
theory. 

It is no wonder, then, that meteorological 
theory in many cases has degenerated into 
pretty pieces of mathematical exhibition, where 
the postulates lack all resemblance to the con- 
ditions actually found in the air and where the 
results cannot be checked, because they belong 
to an imaginary world where all equations are 
linear, the atmosphere homogeneous and _in- 
compressible and the surface of the earth as flat 
as in pre-Ptolemeian days. This flight of meteoro- 
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Comments on Meteorological Research 


C.-G. Rosssy, Massachusetts Institute of Technology 
(Received November 1, 1933) 


logical theory from reality has caused such a 
deep distrust, particularly in the United States, 
of theoretical investigations that synoptic mete- 
orologists have restricted themselves to an accu- 
mulation of weather-map experience which is 
seldom or never interpreted except in the most 
superficial sense. Education in even elementary 
classical physics is not considered necessary or 
even desirable. This development has gone so 
far that even well-established, simple physical 
principles, such as dynamic heating or cooling by 
vertical displacements, are often completely dis- 
regarded in the analysis of synoptic develop- 
ments. 

It would seem, then, that the principal task 
of any meteorological institution of education 
and research must be to bridge the gap between 
the mathematician and the practical man, that 
is, to make the weather man realize the value of 
a modest theoretical education and to induce the 
theoretical man to take an occasional glance at 
the weather map. 

At this point it is important to emphasize that 
the purpose of meteorological theory for many 
years yet will be far more limited than that of a 
theory in, say, celestial mechanics, the object of 
which is to predict future events. In meteorology 
we are satisfied if the theories help us to under- 
stand the individual atmospheric phenomena 
which combine to make the weather. Numerical 
prediction on the basis of a reasonably complete 
physical theory is too far off to deserve immediate 
attention. 

Because of this modest goal, the restrictions 
and postulates of theoretical meteorology often 
constitute a first step which is both permissible 
and desirable, provided they do not affect the 
essential character of the problem. However, 
such an intelligent adaptation of the problem 
is possible only through close cooperation be- 
tween the mathematician and the synoptic 
meteorologist. 
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The polar front theory, developed by V. and 
J. Bjerknes, Solberg, and Bergeron, beyond a 
doubt represents the most successful effort as 
yet made to bridge the gulf that sepgrates the 
two meteorological camps. 

The introduction of air masses and fronts, 
conceptions based on the occurrence of large air 
currents of horizontally approximately homo- 
geneous characteristics, separated by narrow 
transition zones of steep gradients and strong 
atmospheric activity, has given to the weather 
map a living reality which it entirely lacked in 
the days of moving isobars, highs and lows. 

Although the experienced man often may 
recognize the fronts by a glance at the map this 
is not always the case. Often they are masked by 
friction and radiation influences near the ground. 
“False” stationary fronts may appear due to 
abrupt changes in the characteristics of the 
underlying surface (for instance, at a coast line). 
But even if we allow for the uncertainty or errors 
introduced due to the fact that our surface data 
are not always representative of conditions in 
the free atmosphere, there is such a difference 
between the analyses made in the early part of 
the last decade with those of the present day 
that one is led to the conclusion that the polar 
front theory is just as much a method of looking 
at the weather map as an actual statement of 
physical laws. For this reason it would be absurd 
to claim that the polar front theory does not 
apply west of the Rocky Mountains just because 
ideal cyclones are not generally found on the 
weather maps for this region. They cannot be 
found, partly because the network of stations is 
entirely too wide-meshed in the western part of 
the United States and partly because the fronts 
are rather diffuse in the surface layers on account 
of the almost identical oceanic trajectories of 
most of the participating air currents. The polar 
front theory does not state that ideal cyclones 
should be found along the Pacific Coast. It does 
say, however, that it is always possible to de- 
scribe the structure and the behavior of extra- 
tropical cyclones in terms of air masses and fronts. 
Analyses made on airway meteorological ma- 
terial from the San Joaquin Valley and Southern 
California, a region which is far more sluggish 
meteorologically than the state of Washington or 


British Columbia, show that this is entirely 
feasible. 

In the theoretical investigations by the Nor- 
wegian School the narrow transition zones have 
been supplanted by surfaces of discontinuity. 
Thus the theoretical problem resolves itself into 
a study of the dynamics of internal atmospheric 
boundaries. V. Bjerknes and Solberg have studied 
various types of wave motion in such boundaries. 
Though the restrictions introduced to facilitate 
the mathematical treatment are too great to 
permit an immediate acceptance of their findings 
as characteristic of the real atmosphere, they 
have obtained several results of definite practical 
significance, such as values for the rate of 
progress of these waves and a criterion for the 
stability of front surfaces of varying angles of 
incidence. 

The existence of rapidly progressing waves of 
dimensions similar to those of extra-tropical 
cyclones was first demonstrated synoptically by 
Bergeron in 1924 as a result of a painstaking 
study of a certain European synoptic situation, 
with the aid of an enormous observation ma- 
terial. With the improved technique of front and 
air mass analysis it has become a simple matter 
to identify some of these waves. The daily 
working maps on file at the Massachusetts 
Institute of Technology contain a number of 
cases of quasi-stationary fronts with progressive 
waves, although these maps have been analyzed 
in much less detail on account of the limited 
time available. 

The vertical structure of an air mass depends 
partly upon its source region, partly upon its 
trajectory. If these two elements are known, it 
is possible to predict, in a general way, the 
vertical distribution of temperature and hu- 
midity. It is possible to check these conclusions 
in a qualitative fashion through a study of the 
typical cloud forms associated with each air mass, 
since each type of temperature and humidity 
stratification has its characteristic cloud forms. 
A detailed check, however, is possible only from 
a study of pressure, temperature and humidity 
records obtained through airplane or kite sound- 
ings. Using a large amount of kite data from the 
Weather Bureau aerological stations my col- 
league, Dr. Willett, has recently completed an 
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investigation of the vertical structure of Amer- 
ican air masses. The results of this work will be 
given in the next issue of this Journal. A pre- 
liminary study of the air masses of the Pacific 
based on airplane records from Honolulu, San 
Diego and Seattle, has been made at M.I.T. by 
Byers and will soon be published by the Scripps 
Institution in La Jolla. However, the amount of 
aerological material from the West Coast and 
from the North Pacific is too scanty for con- 
clusive discussion. 

Parallel with these synoptic investigations, the- 
oretical studies have been made which have led 
to the development of a simple graphical method 
of identifying air masses with the aid of their 
most conservative properties, entropy of the 
dry air and specific humidity. An atmospheric 
sounding, if plotted with these two quantities as 
coordinates, gives a characteristic curve which, 
at least during the winter season, remains re- 
markably constant throughout the life history 
of the air mass. 

Outstanding problems still awaiting their final 
solution in connection with the study of air mass 
properties are the following: 

1. A determination of the rate at which an 
air mass from a given source region picks up 
moisture from the ground and warms up because 
of increased radiation heat income, coupled with 
increased turbulent or convective transport. The 
theoretical solution of this problem should be 
checked through a study of aerological ascents 
along the trajectory of the air mass. 

2. A study of the rate of accumulation of 
energy available for convection and thunder- 
storms in stagnant air masses, coupled with a 
study of the mechanism through which this 
energy is released. 

Both problems are of significance in con- 
nection with summertime forecasting and both 
require a large amount of aerological material 
for their solution. Such material is not now 
available and would require the organization of a 
small-mesh aerological network, each station 
making several soundings a day. 

An aerological airplane station was established 
at Boston in November, 1931 and daily flights 
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are made up to a level of approximately 5 km. 
The data from this station, combined with the 
soundings from the Weather Bureau aerological 
stations in Cleveland, Chicago and Omaha, all 
at about the same latitude and fairly evenly 
spaced, offer an excellent opportunity for a 
charting of the vertical structure of passing warm 
sectors and polar air wedges. A preliminary 
survey of the available observation material 
seems to indicate that in general our cold air 
masses are more shallow than those of Europe, 
a result which becomes understandable when one 
considers the great difference in latitude. Thus 
Scandinavian meteorologists, in discussing Euro- 
pean conditions, sometimes dismiss the first 
kilometer as non-representative of the air mass, 
while in the United States, especially in the 
South and Southeast, a stagnant polar air mass 
sometimes does not reach higher than fifteen 
hundred meters. 

During the last few years attention has been 
focussed upon the strong pressure and tempera- 
ture changes at the base of the stratosphere. 
Bjerknes and Palmén interpret these variations 
as wholly or in part dynamic effects of the 
succession of cold and warm currents within the 
troposphere. On the other hand, Austrian mete- 
orologists have claimed that the substratosphere 
must be regarded as an additional center and 
source of atmospheric activity. If the Scandi- 
navian notion is correct, one would expect that 
the stratospheric variations in temperature and 
pressure over the United States would be less 
marked because of the fact that our wedges of 
polar air do not reach so high up. 

The preceding remarks should suffice to dem- 
onstrate the fecundity of the Norwegian ideas as 
a working hypothesis. American participation in 
the research work which has sprung up with the 
polar front theory is still very slender and 
largely confined to a sterile discussion of the 
somewhat antiquated Norwegian conceptions of 
the early ‘twenties. Vigorous American participa- 
tion in research on current problems is desirable 
but must be preceded by general acceptance of 
the modern methods of front and air mass 
analysis. 
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Fluid Flow Near a Smooth Surface 
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T is known from model experiments that the 
air in contact with the surface of a streamline 
body moves with it, and that the regions of easy 
flow which surround the body are connected to 
the surface through a relatively thin layer of air, 
within which the velocity falls to zero. It is 
within this boundary layer that the resistance to 
motion is transmitted to the surface. 

At speeds of flight, the flow in the boundary 
layer over the front part of a streamline body is 
regular, or as it is more often called, laminar, 
and the elements of fluid passing through any 
point in the layer behave in a simple and definite 
manner. At some distance beyond the nose of 
the body this regular motion breaks down, the 
movements of the elements become confused, and 
the flow is said to be turbulent. 

Measurements made with small Pitot tubes 
and hot-wire meters have shown that the velocity 
in a turbulent boundary layer falls rapidly as 
the distance from the surface diminishes. Vis- 
cous forces then predominate, and it has been 
suggested that the flow simulates more and more 
closely a streamline character as the surface is 
approached. The intensity of surface friction 
would then be given by the product of the 
velocity gradient at the surface and the coeffi- 
cient of viscosity. To determine whether this 
relation held, Stanton! made very careful meas- 
urements at the surface of a pipe in which the 
flow was turbulent. It was known that a reliable 
measurement of the velocity gradient at the 
surface could not be obtained with a Pitot tube 
of the usual type, even when the diameter of 
the mouth was very small, since a sufficiently 
close approach to the surface would not be 
possible: or with an exceedingly fine hot wire, 
because in addition to the forced convection due 
to the wind stream there would be a loss of 
heat by conduction across the thin layer of air 


1 Stanton, Marshall, Bryant, Proc. Roy. Soc. A97 (1920). 
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between the hot wire and the surface. To over- 
come these difficulties, an ingenious surface tube 
(Fig. 1), which had for its inner wall the surface 
of the pipe itself, was designed and used by 
Stanton. Because of the extreme smallness of 
the tube, a special method of calibration had to 
be devised to determine the effective distance 
from the surface corresponding to the speed 
calculated from the pressure measured at the 
mouth. The product of the measured velocity 
gradient at the surface and the coefficient of 
viscosity of the air was found to be in fairly 
close agreement with the intensity of surface 
friction. More recently, the intensity of friction 
on the surface of an aerofoil*® has been de- 
termined at the National Physical Laboratory 
from measurements of velocity gradient taken 
with a particular form of Stanton tube con- 
structed on the top of a circular rod designed to 
pass with a small clearance through holes drilled 
in the polished surface of a metal model. Three 
tubes were used. Photographs of the mouths of 
these tubes are given in Fig. 1. An interesting 
feature shown by the calibrations of these tubes 
is that the ratio of the effective distance to the 
width of opening increases as the width de- 
creases. The effective distance of the smallest 
tube (No. 1) is in fact greater than the width 
of the opening (0.0020 inch), and the observa- 
tions taken with this tube were not appreciably 
nearer the surface than those taken with the 
largest tube No. 3 (0.0044 inch). The distribution 
of frictional intensity was determined from meas- 
urements of velocity taken at 0.002 to 0.003 
inch from the surface. The frictional drag of the 
aerofoil obtained from an integration of the 
frictional intensities was found to be in fairly 


* The intensity of the surface friction on a circular 
cylinder has also been measured in a similar manner. 
(R. and M. No. 1369.) 

? Fage and Falkner, Proc. Roy. Soc. A129 (1930). 
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close agreement with predictions by two other 
methods. 

The work with Stanton tubes leads then to 
the conclusion that in eddying flow the gradient 
of mean speed tends to a value at the surface 
which seems to be the same as that for a stream- 
line motion having the same intensity of friction 
on the surface. Is then the flow in a surface 
layer of small but finite thickness of a streamline 
character? To obtain an answer it is necessary to 
watch the movements of particles very near the 
surface. 

The scale of turbulence near a surface is 
small, and it is most important that the method 
of observation should not interfere in any manner 
with the internal motions of the fluid. The 
method of examining fluid motion with an ultra- 
microscope*® which has been in use for some time 


’ Fage and Townend (with Appendix by Taylor), Proc. 
Roy. Soc A130, 656 (1932). 
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at the National Physical Laboratory satisfies this 
condition, since neither extraneous particles nor 
measuring instruments are introduced into the 
fluid. The principle of the ultramicroscope de- 
pends on the fact that very small particles always 
present in fluids, such as air or tap water, but 
invisible in ordinary light even under the most 
powerful microscope, become visible when in- 
tensely illuminated and observed against a dark 
background. Even particles smaller than the 
wave-length of light become visible as bright 
points of light, if the intensity of the light beam 
is sufficiently great. Fig. 2 gives a photograph of 
the ultramicroscope arranged to examine tur- 
bulent flow in a square pipe, with 0.89 inch 
sides, a flow which is of the same kind as 
that in the turbulent boundary layer of a flat 
plate or a streamline body. A beam of intense 
bright light passing through a window in one 
side of the ; .pe illuminated particles contained 
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in the water, and the motions of these par- 
ticles as they flowed downstream were viewed 
through an adjacent window with a microscope. 
When the motion was turbulent the particles 
passing through any fixed point moved in 
different directions, and the short traces of the 
paths seen in the microscope appeared as bright 
rectilinear streaks inclined at various angles to 
the mean direction of flow. The maximum 
obliquities of these streaks were measured. Meas- 
urements were also made of the mean speed and 
the maximum velocity fluctuations in the axial 
direction, by a technique which was based on the 
fact that illuminated particles were brought to 
apparent rest by motion of the objective of the 
microscope at the appropriate speed. 

It was found that the maximum obliquity of 
particles in an axial plane normal to a wall of 
the pipe increased slowly with the distance from 
the axis until a maximum value was reached about 
0.05 inch from the wall. The value then fell steadily 
to zero at the wall, so that particles seen when 
looking along and very close to the wall appeared 
to be moving in laminae parallel! to the wall. The 
maximum obliquity in planes parallel to the wall 
increased progressively as the wall was approached 


and at the small distance of about 0.003 inch the 
large angle of 70° was observed. On viewing the 
surface with a microscope having a very high 
magnification (200) illuminated particles within 
a surface layer of thickness 0.001 inch were 
brought into focus. The faster moving particles 
in this layer were obviously moving in sinuous 
paths, and careful observation showed that the 
slower ones, some of which were estimated to be 
not more than 2.5X10-* inch from the surface, 
crossed and recrossed several times a hairline in 
the eyepiece of the microscope. Occasionally, 
when a group of the faster particles moved 
laterally sideways the slower particles would also 
shift laterally. Also, the axial motions of the 
slower particles were often jerky. 

Fairly conclusive evidence was therefore ob- 
tained to show that when the general motion is 
turbulent the flow in a very thin surface layer is 
not regular, although the friction on the surface 
can be determined fairly closely from the relation 
f=u(dV/dy),-0, where V is the mean velocity at 
a distance y from the surface, and the value of 
(0V/dy),-0 is predicted from values of V meas- 
ured with a small surface tube. 

The ultramicroscope has also been used to 
observe the nature of the eddying wake behind 
long cylindrical bodies.* This work leads to the 
interesting conclusion that the disturbances in 
the eddying wake behind any long cylindrical 
body, bluff or streamline, placed with its length 
normal to the direction of flow, will be three- 
dimensional, under conditions of flight, even 
although the mean flow in planes at right angles 
to the length tends to be two dimensional. It is 
perhaps of interest to mention that photographs 
have recently been taken of some interesting 
views of fluid flow revealed with an _ ultra- 
microscope. 


4 Fage, A.R.C. (R. and M. No. 1510). 
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Stresses in Wire Wheels 


A. J. Sutton Pipparp, Bristol, England 
(Received October 10, 1933) 


URING the last two or three years the 
writer has been carrying out a research into 
the stresses caused in a wire wheel under loads 
applied to the rim. This work has been done for 
the Aeronautical Research Committee and the 
results have been published in a series of papers 
which are enumerated at the end of this note 
which gives a brief resumé of the problem, 
method of attack and certain results. 

A wire wheel consists of a comparatively rigid 
hub and a rim connected to it by a large number 
of small section spokes which may be either 
radial or inclined to the radius. Since the flexural 


rigidity of the spokes is small it is sufficiently ’ 


accurate for analytical purposes to assume that 
they are pin jointed to the hub and rim. With 
this assumption it is, of course, impossible for the 
radially spoked wheel to resist actions on the rim 
which are tangential, e.g., those due to the appli- 
cation of brakes: for efficient resistance to such 
forces the spokes must be inclined in the plane 
of the rim. To resist side forces on the rim the 
spokes must also have an inclination in a plane 
perpendicular to that of the rim. 

The problem to be considered, therefore, was 
that of a wheel having radial or inclined spokes, 
held at the hub and subjected to loads at the 
rim acting radially, tangentially and trans- 
versely. An attack was first made by means of a 
strain energy treatment but the resulting equa- 
tions were too complex for a complete solution 
to be obtained. This analysis did, however, yield 
certain important relationships between the 
resultant actions which were required. A second 
line of attack was suggested to the author by 
Professor R. V. Southwell, which was instru- 
mental in completing the solution of the general 
problem. 

It was assumed that the spokes could be 
replaced by a continuous disk which had the 
peculiar property of transmitting load in one 
direction only. If the spokes were radial, the disk 
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could have stresses only in a radial direction; if 
they were inclined, the disk transmitted stresses 
only along the inclined lines. If more than 
one system of spokes were used in the wheel, 
as is commonly the case, each system of spokes 
was assumed to be replaced by an appropriate 
disk. The displacements of the rim under load, 
at any point y from an axis of reference are 
denoted by wu and 7 in the radial and tangential 
directions, respectively. R is the radius to the 
centroid of the rim, Ro, the radius to the tangent 
circle on the hub, J, the moment of inertia of the 
cross section of the rim about its axis of bending, 
a, the cross-sectional area of the rim, F,, the 
modulus of elasticity of the rim material, E2, the 
modulus of elasticity of the disk material, ¢,, the 
inward pull per unit length of rim of the disk 
acting at ¢ to the radius, fs, the inward pull per 
unit length of rim of the disk acting at — « to the 
radius, where « is the angle of inclination of the 
spokes to the radii in the plane of the wheel. If 
the thickness of each disk at the rim is A, we have 


t, = (E£,A/R)(u—v tan e); 
to = tan e). 


Also let 7 be the tangential force on the rim at y, 
V thé radial force and M the bending moment. 
Then 

H = —(E,a/R)(u+dv/dp) (1) 


M = (2) 
The differential equation for rim displacements is 
then 
2d'u/dyi+(1+K, cos €)(0?u/dy) 
— Kyu sin e tan e=0, 


where K, The complete solution 
of this equation is 
u=A sinh y~¥+B cos ay sinh By 
+Csin ay cosh BY+D cosh yw 
+E cos ay cosh BY+F sin ay sinh By, 
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y= (X-2/3)4, 


1 Ky, sece 
M=—+———- (2+sin’? e), 
27 6 


where 


Ki K 2 3 
+K,{2+sin? (27/4) tan? e)}+sec ; 


Since for inextensional strain of the rim 
u= —du/dy, 


a complete solution of the displacements of the 
rim is found and hence the forces ¢; and fy at 
any point are known. In the limit when ¢«=0, 
y is zero, and a and B reduce to the values for the 
radially spoked wheel, viz.: 


The constants A, B, C, D, E and F are evaluated 
from the known conditions of the problem, and 
the work, although somewhat lengthy, is not 
difficult. Having found ¢; and fe for the disk, the 
load in a spoke of the actual wheel is determined 
approximately by integrating these stresses over 
a segment of the disk replacing the spoke con- 
sidered. It will be apparent that the greater the 
number of spokes in the actual wheel the better 
will be the approximation. 

In order to check the accuracy of the solutions 
thus obtained, a series of experimental wheels 
were constructed. The rims, four in number, of 
different sizes, were cut out of mild steel plate 
and the spokes were made of 4 B.A. aeroplane 
tie rods and 3” rods. Each rim was used with 
three, four, six, twelve and twenty-four spokes of 
each size and the loads in the spokes were meas- 
ured by accurate extensometers for comparison 
with the calculated values. From these measure- 
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ments the values of the bending moment in the 
rim could be deduced from the strain energy 
analysis and could thus be compared with the 
more approximate solution. 

The general agreement between calculated and 
measured values of loads and rim actions was 
excellent and the experiments showed con- 
clusively that the approximate solutions obtained 
were quite satisfactory for any wheel with more 
than 30 spokes. The formulae deduced for spoke 
loads, etc., were rather complicated and in view 
of the good experimental results obtained, these 
formulae were further simplified for convenience 
of calculation with very little loss of accuracy. 
The simplified formulae are as follows: Let 


P be the inward radial load on the rim, 
Ty, the maximum compression in any spoke, 
Ty, the maximum tension in any spoke, 
mR, the maximum bending moment in rim. 


Then Ty/P=1—e-*?; 
mo’ /P =Bo/2(Ki)}. 


These results are only applicable to the case of 
the wheel with radial spokes. When the spokes 
are inclined the results are too complicated for 
inclusion in this note but in the original papers 
cited below full particulars, together with curves 
for rapid estimation, are given. 

In dealing with the case of transverse loads on 
the rim certain assumptions were made to reduce 
the problem to manageable proportions and it 
was found that by a simple transformation of 
some of the quantities involved the stresses in 
the spokes and actions in the rim could be deter- 
mined from curves already plotted for the radial 
and tangential cases of loading. 
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Effects of Compressibility on High Speed Flight 


Joun Stack, National Advisory Committee for Aeronautics, Langley Field, Virginia 
(Received October 28, 1933) 


T has been known that the air forces exerted 

on aerodynamic bodies follow the so-called 
“square law”’ except for the effects of scale and 
compressibility. The effects of dynamic scale, 
or Reynolds number, have become rather well 
known during the last 10 years, largely because 
of research conducted in the N.A.C.A. variable- 
density wind tunnel. The effects of compressi- 
bility have commonly been neglected because 
until the relatively recent development of the 
last Schneider trophy aircraft the speeds have 
been low as compared with the velocity of sound, 
and the consequent local pressures over the 
surfaces of high speed airplanes have differed 
but slightly from atmospheric pressure. At the 
present time, however, the speeds associated with 
the fastest airplanes approach 60 percent of the 
velocity of sound, and the induced velocities 
over their exposed surfaces lead to local pressures 
that differ appreciably from the pressure of the 
atmosphere. When this condition exists air can 
no longer be regarded as an incompressible 
medium. The effects of compressibility on the 
aerodynamic characteristics of airfoils have been 
under investigation by the N.A.C.A. in the high 
speed wind tunnel, and it is the purpose of this 
paper to examine the possibility of further 
increases in speeds in the light of this relatively 
recent research. 

The following computations are made for a 
hypothetical airplane which, however, is not 
beyond the limits of possibility. The airplane is 
a scaled-up model which has been tested in the 
variable-density wind tunnel as part of an 
extensive investigation of wing-fuselage inter- 
ference, and represents one of the best wing- 
fuselage combinations thus far produced. The 
geometric scale factor applied to the model 
dimensions was determined from the ratio of the 
model fuselage diameter to the diameter of a 
fuselage sufficiently large to house a 2300- 
horsepower Rolls-Royce R type engine. The 


model from which the airplane is derived is a 
mid-wing cantilever monoplane having a fuselage 
of circular cross section, and a straight tapered 
wing having the N.A.C.A. 0018 section at the 
center and the N.A.C.A. 0009 section at the tip. 
The principal full-scale airplane dimensions are 
given in Table I. An outline drawing of the 
airplane for which computations will be made is 
given in Fig. 1. 


TABLE I. 
Wing chord (average) . 4.85 ft. 
i 6 


Aspect ratio 
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Fic. 1. Hypothetical high speed airplane. 


The model, without tail surfaces, tested in the 
variable-density wind tunnel at a Reynolds 
number of 3,000,000, was found to have a 
minimum drag coefficient of 0.0115. If 10 percent 
of this value be added to account for tail-surface 
drag the resulting coefficient is 0.0127. No drag 
allowances need, or should be, considered for 
landing gear, cooling systems, or other appen- 
dages. The airplane should have completely 
retractable landing gear or be catapulted and 
landed on water with the aid of suitable retract- 
zble hydrodynamic surfaces. The surfaces of the 
airplane must be free from protuberances or 
excrescences because of the large amount of 
power required at very high speeds to overcome 
the drag of even the slightest surface irregularity. 
The engine will be completely enclosed and 
cooling can be provided by skin-type radiators. 
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Vision for the pilot can be obtained by trans- 
parent fuselage sections or indirectly by mirror 
systems. Inasmuch as the condition of the air- 
plane herein assumed is not impossible of attain- 
ment, the airplane may be considered as an 
example of a high degree of aerodynamic refine- 
ment but not as an “ideal” airplane. 

If the effects of compressibility are, for the 
moment, neglected, the computation of the 
maximum speed when a 2300-horsepower engine 
is used becomes relatively simple. The drag 
coefficient can be considered as less than the 
variable-density wind-tunnel value because of 
the higher Reynolds number of the airplane. 
The minimum profile drag coefficient may be 
reasonably estimated as 


R.N. 
R.N. flight 


This assumed reduction in profile drag is based 
on the results of tests on airship and airfoil 
models over a wide scale range, and is probably 
not greatly in error. In order to estimate the 
full-scale drag coefficient a speed of 550 miles 
per hour is assumed and the corresponding 
Reynolds number is computed. The full-scale 
minimum profile-drag coefficient then becomes, 
according to the preceding formula, 0.0098. A 
reasonable assumption regarding the probable 
propeller efficiency must next be made. The 
propeller pitch will exceed the optimum value, 
but the decrease of efficiency from this source 
will probably not be very great. The tip speed 
will be high but losses from this source may be 
kept to a minimum by suitable gearing and a 
reduced propeller diameter. The large forward 
thrust of the propeller will necessitate wide 
blades and probably a multibladed propeller. 
There are, of course, no direct experimental data 
on which to base an accurate estimate of the 
efficiency for such propellers, but it is believed 
that efficiencies at least as great as 75 percent 
may be obtained. 

Assuming, then, 75 percent for propeller 
efficiency the maximum speed of the airplane, 
neglecting the induced drag, becomes 


Cp, min = Cp, min(tunnel) ( 


Vinax = (2 hp. 7550/Cp, min ; A)'=830 ft. /sec. 


or 


566 miles pei hour. 


The slight error in the assumed Reynolds number 
from which Cp, min was computed is negligible 
and within the probable error of the assumed 
effect of scale. The induced drag will cause but 
slight reduction in speed. Assuming 2.5 pounds 
per horsepower, the weight of the airplane is 
5750 pounds and the lift coefficient at which 
the high speed flight occurs is approximately 
0.05. The induced-drag coefficient, Cp,, is there- 
fore 0.0001, approximately 1 percent of the total 
drag. 

It is interesting, now, to compute the landing 
speed of the airplane. The maximum lift coeffi- 
cient as determined by tests in the variable- 
density wind tunnel is 1.516. The landing speed 
computed from this value is therefore 


V=(2L/pC,A)!=150.6 ft./sec. 
or 
V =103 miles per hour. 


The effects of compressibility as indicated from 
airfoil tests in the high speed tunnel will ma- 
terially affect this calculated performance. Re- 
cent tests' demonstrate that the effects of com- 
pressibility are of considerable importance for 
speeds above 40 percent of the velocity of sound 
(306 m.p.h.). This result has been corroborated 
by tests on a series of symmetrical airfoils, the 
results of which are now being prepared for 
publication. The variation of the minimum drag 
coefficient with speed for a symmetrical airfoil, 
the N.A.C.A. 0012, is shown in Fig. 2. The drag 
coefficient is given as D/}pV?A. The Reynolds 
number for these tests is, however, much lower 
than that for the airplane under discussion 
because of the small size of the model. Unfortu- 
nately, the effects of large Reynolds number 
differences at these high speeds are unknown. 
These data, however, represent the best available 
information and will be applied without scale 
correction. 

The landing speed of the airplane may be 
expected to be unaffected by compressibility 
because at such low speeds (103 m.p.h.) no 
measurable compressibility effect on maximum 
lift has been observed. It may therefore be 


‘John Stack, The N.A.C.A. High-Speed Wind Tunnel 
and Tests of Six Propeller Sections, Tech. Report No. 463, 
N.A.C.A., 1933. 


42 JOHN STACK 


-10 


S 


| 
| 


| 


Minimum Profile —Drag Coefiicient, Co, otis 


| Velocity of Sound 
| 

0 | 

0 200 400 600 800 

Speed, Miles per Hour 


Fic. 2. Effect of compressibility on minimum profile-drag 
coefficient. N.A.C.A. 0012 airfoil. 


taken that the landing speed of the airplane is 
not beyond reason. 

The high speed of the airplane will, however, 
be materially reduced because of the effects of 
air compressibility. These effects have been esti- 
mated from the results of tests of the N.A.C.A. 
0006, 0009, and 0012 airfoils in the high speed 
wind tunnel. In this estimation the drag co- 
efficient of the airplane, 0.0127, was corrected 
by the previously given formula for scale effect 
to a Reynolds number corresponding to 306 
m.p.h. Above this speed the drag coefficient was 
corrected in accordance with the high speed 
tunnel data. The resulting minimum-drag curve 
for the airplane is shown in Fig. 3. 

A direct computation of the maximum speed 
cannot be made because of the Cp, variation 
with speed. An estimate could be made by a 
series of successive approximations but a better 
method which shows the compressibility effects 
in terms of horsepower differences for a given 
speed or in terms of speed for a given horsepower 
is to compute the horsepower-required curves. 
In this computation the induced drag and the 
additional profile drag due to lift have been 
added to the minimum profile drag of the wing. 
The weight has been held constant at 5750 
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pounds. The power curves are shown plotted in 
Fig. 4. 

Assuming, as before, a propeller efficiency of 
75 percent and 2300 horsepower, the speed loss 
due to the compressibility effect is shown by 
curves A and B as 39 m.p.h. Expressed in terms 
of horsepower at 524 m.p.h., the maximum speed 
corrected for compressibility, the effect is 325 
horsepower, sufficient power to propel the air- 
plane at 300 m.p.h. 

The increasing importance of the compressi- 
bility effects has been foreseen by several re- 
search agencies, but the largest amount of work 
has been done by the National Advisory Com- 
mittee for Aeronautics. Early tests? indicated 
the trend of aerodynamic characteristics with 
speed and were followed by improved means for 
studying the phenomena. 

Some results are already available to designers 
in the form of new wing sections developed in 
the high speed wind tunnel. The drag of one of 
the new sections, the N.A.C.A. 204, is con- 
siderably lower than the drag for the section of 
our tapered wing having the same thickness but, 


2L. J. Briggs, G. F. Hull and H. L. Dryden, Aero- 
dynamic Characteristics of Airfoils at High Speeds, Tech. 
Report No. 207, N.A.C.A., 1925. 
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most important for purposes of high speed flight, 
the speed at which the compressibility burble of 
the N.A.C.A. 204 occurs (evidenced by the rapid 
rise of the drag coefficient at high speeds) is 
delayed to higher speeds. Unfortunately, no data 
are as yet available on the maximum lift and 
the effects of thickness changes for the N.A.C.A. 
204. However, making reasonable assumptions 
regarding the effects of thickness variations and 
applying the indicated drag reduction to the 
parts of the airplane affected by the change to 
the new wing section, the high speed portion of 
a new power required has been computed and 
plotted as curve C in Fig. 4. 


A comparison of curve C with curves A and B 
indicates a decided improvement at the lower 
speeds, due, of course, to the lower drag of the 
new section. A considerable increase in the 
maximum speed is also indicated. For the as- 
sumed power and propeller efficiency the maxi- 
mum speed of the airplane is increased from 524 
to 544 m.p.h. In terms of horsepower at 544 
m.p.h., a direct saving of 315 horsepower is 
indicated. The improvement at the extremely 
high speed can be considered as due to the 
delayed compressibility burble. 

It is apparent from the foregoing computations 
that much higher speeds than those so far 
attained are possible and likely. The attainment 
of these higher speeds is, however, dependent to 
a large extent on the effects of compressibility, 
first because of the compressibility effect on the 
airplane drag, and second because of the com- 
pressibility effect on the propeller efficiency. 
The aerodynamic characteristics of propellers 
suitable for such high speeds must also be 
investigated. The range of pitches for which tests 
have already been made must be extended and 
the effects of wide variations in blade shape and 
number of blades must also be studied. 

The study of the compressibility phenomena 
is of relatively recent origin and much remains 
to be done. Practically nothing is known of the 
combined effects of scale and compressibility at 
very high speeds, nor is anything known of 
aerodynamic interference at very high speeds. 
The success so far achieved, however, is en- 
couraging and it is reasonable to expect further 
aerodynamic improvement as the knowledge and 
understanding of the nature of the compressi- 
bility phenomena expand. 
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Speed of Air Driven Rotors in Gyroscopic Instruments 


A. F. Sprtuaus, Massachusetts Institute of Technology 
(Received October 10, 1933) 


HE importance of blind flying instruments 
is now generally accepted and gyroscopic 
instruments are used very widely in this con- 
nection. Although European design has employed 
electrically driven as well as air driven rotors, 
it is almost without exception the practice in 
this country to use the air jet drive which, by 
virtue of its simplicity in the instrument itself, 
is very desirable and reliable. The usual arrange- 
ment of such a system is to keep the whole 
inside of the instrument casing at a reduced 
pressure thus causing the differential pressure 
between the atmosphere and the reduced pres- 
sure in the casing to force air through the jet to 
impinge on the rotor buckets. 

The speed at which the rotor is run in these 
instruments is of great importance. Excessive 
rotor speed may prove detrimental to the 
bearings while, if the speed of the rotor drops too 
far below its normal value, the nutations of the 
gyroscope begin to be troublesome and the 
instrument tends to become generally sluggish. 

It is the purpose in the following simple 
considerations to develop an expression which 
will predict the rotor speed at various atmos- 
pheric conditions and thus anticipate the per- 
formance of the instrument at those conditions 
for the particular suction value available. 

The rotor speed will depend on the suction on 
the instrument and also on the conditions of 
temperature and pressure in the atmosphere. 
Then when the rotor is running steadily, 


The Jet Impulse Torque = Windage Resistance Torque 
+Bearing Friction Resistance Torque. 


At the high speeds at which these rotors run the 
bearing frictional torque will be mtch the smaller 
of these two factors, the windage or air resistance 
playing the predominant part. 

Considering first, therefore, the drop in rota- 
tional speed due to windage, if marop is the drop 
in revolutions per sec., 
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Narop =a function (p, n, 7, 4), 


where p is the air density, m is the rotational 
speed, 7 is the effective radius and uy is the air 
viscosity. 

Expressing “drop in a power series, 


Narop = + Ci etc. 


Considering this for dimensional homogeneity, 
the expression 


Narop = Bn X function of (u/ pur?) = Bng(R,), 


where B is a constant, ¢ is some function, and 
pnr?/u= Reynolds number. 

If there were no resistance, the rotational 
speed would be equal to u/2mr, where wu is the 
velocity through the jet. 

The drop in revs., therefore 


u/2xr—n=Bno(R,) 


+drop in revs. due to bearing frict. 


—drop in revs. due to bearing frict. 


Now considering bearing friction, suppose 
there is a constant frictional torque, Fr; then for 
steady rotation (considering bearing friction 


only) 
Copu?r = Fr, 


Copu? representing the effective momentum of the 
jet. This whole term will be smaller than the 
windage term in the complete expression and 
therefore very little error will be introduced by 
writing 
C3pn?r? = F 
n=K,/p} 


(assumption is that wu is proportional as a first 
approximation to ”). Marop due to bearing friction 
may, therefore, be written 


Narop trict. = K/p'. 
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Fic. 2. Sperry Directional Gyro. Rotor speeds at various 


air pressures at two different pressure altitudes. 


The complete expression is then 
n(1i+Bo(R,) 


and when the differential pressure is Ap through 
the nozzle, then u=a(2Ap/p)?. 

In this formula for the rotor speed, therefore, 
only a, K and B¢(R,) have to be evaluated. 

a is simply the orifice coefficient of the nozzle 
through which the air issues to impinge on the 
rotor. @ is a function of Reynolds number but 
its variation is small enough to be neglected; 
a has been assumed here = 0.95. 

K is the constant containing the bearing 
friction effect and may be evaluated in two 
ways. It is evident that when = zero 


K =(a/2znr)(2Ap)}. 


Thus by finding the differential pressure at 
which the rotor stops K can be found. 
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Alternatively, if two curves are available of 
n vs. Ap for two different values of p, then two 
points having the same value of R, can be 
picked. If the conditions at these points be 
mApip; and and pe 


2Api\! K 
py 

a (=) K 

pe n2(p2)* 


This may be solved for K. When K is deter- 
mined, from the curves of m vs. Ap values of 
Bo(R,) can be calculated for various values of 
Re 

Whereas K will vary from instrument to 
instrument of the same type because of slight 
variations in the bearings, it is evident that 
B¢(R,,) will be substantially the same for all 
instruments of one type. Figs. 1 and 2 show the 
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Sperry Directional Gyro 
Fic. 4. 


variation of rotor speed with suction for various 
conditions of pressure and temperature for the 
Sperry Artificial Horizon and Directional Gyro, 
respectively. The values of Bd¢(R,) calculated 
from these curves are shown plotted against 
(R,,) in Fig. 3. 

Examination of Fig. 3 shows that the windage 
factor may be expressed over the practical range 
as B(R,,)* where 6 is the gradient of the line on 
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the logarithmic plot; also that the windage factor 
for the Directional Gyro is much greater than 
that of the Horizon and that 6 is small for the 
former instrument. If 6 were zero, then it is 
evident that windage loss would be entirely turbu- 
lent; thus the value of 6 indicates the degree to 
which viscosity is entering into the resistance 
effect. 

5is, therefore, greater for the Artificial Horizon 
because the rotor is completely enclosed in a 
casing with small clearances. This casing also 
has the effect of reducing the total value of 
Byor.(Rn)* for the Horizon as it reduces the size 
of the turbulent eddies. In the Directional Gyro 
the rotor is open and the turbulent eddies are 
large, viscosity enters less, and, therefore, 6 is 
small and the value of Bpir.(R»)> is much greater. 
This indicates well the beneficial effect of en- 
casing the rotor in such an instrument. 

The above analysis permits the separation of 
the various factors affecting the rotor speed and 
it then becomes evident that with an increase 
of altitude with a constant suction on the 
instrument, the decreased pressure will cause the 
rotor to speed up while the decrease in tempera- 
ture may cause K to increase and thus tend to 
slow the rotor down; but with carefully designed 
bearings the net result is an increase of rotor 
speed with altitude. 
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Institute of the Aeronautical Sciences 


AIMs 


HE Institute of the Aeronautical Sciences 

was organized to promote the application 
of science in the development of aircraft. By 
grading members according to their qualifica- 
tions, recognition is given to outstanding pro- 
fessional work and an incentive is provided 
younger members for the investigation of new 
aeronautical problems. The presentation and 
discussion of technical papers at Institute meet- 
ings afford members opportunities to hear and 
meet the leaders in specialized fields of the aero- 
nautical sciences. Through the publication of the 
Journal of the Aeronautical Sciences, members 
are provided with an opportunity to publish and 
to read new knowledge and new applications in 
the field of aeronautical endeavor. By including 
in its membership the leading aeronautical 
specialists in other countries, an interchange of 
international aeronautical scientific progress is 
made possible. The organizers of the Institute 
believe that the future development of aircraft 
depends on scientific research and engineering 
experiment applied to design. As practically 
every scientist and engineer already belongs to 
his own special professional society, in order to 
bring together all specialists contributing to 
aeronautical progress, the dues of the Institute 
have been made nominal. Through the foregoing 
aims, as well as others too obvious to specify, the 
Institute is providing the United States with a 
scientific society similar to those which have had 
such a beneficial influence on aeronautical prog- 
ress in other countries. 


HIsTORY 


The Institute of the Aeronautical Sciences was 
incorporated under the membership corporation 
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law of New York State on October 15, 1932. 
The purpose of the Institute as stated in the 
Charter is as follows: 

“To advance the art and science of aero- 
nautics; to publish works of literature, science 
and art for such purpose; to do all and every act 
necessary, suitable and proper for the accom- 
plishment of any of the purposes or the attain- 
ment of any of the objects or for the furtherance 
of any of the powers hereinbefore set forth, either 
alone or in association with other corporations, 
firms or individuals and to do every other act, 
thing or things, incidental or appurtenant to or 
growing out of or connected with the aforesaid 
science or art, Or powers or any part or parts 
thereof, provided the same be not inconsistent 
with the laws under which this corporation is 
organized or prohibited by the laws of the State 
of New York.”’ 

The incorporators at the first meeting, on 
October 17th, 1932, adopted a constitution for 
the Institute which gives broad powers to the 
Council until such times as the members of the 
Institute wish to make changes. On October 24th, 
the Council adopted By-Laws and elected tem- 
porary officers to serve until the first annual 
meeting of the Institute which was held on 
January 26th, 1933. 

At the Founders’ Meeting held at Columbia 
University on January 26th, 1933, the Constitu- 
tion and By-Laws were approved by the members 
and the temporary Council was re-elected. The 
papers presented at the Founders’ Meeting were 
published and bound copies were distributed to 
all members and to aero clubs, universities, 
libraries and aviation officials in many countries. 

A members’ clubroom and an office have been 
provided, by generous donors, at Rockefeller 
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Center, New York, for the social and adminis- 
trative activities of the Institute. A library of 
books, magazines and pamphlets is now available 
to members. 

The need for such an aeronautical organization 
is best indicated by the membership list which 
includes leaders of aeronautical thought in the 
United States and abroad. 


PUBLICATIONS 


The Journal of the Aeronautical Sciences is 
the official publication of the Institute. It is 
published in collaboration with the American 
Institute of Physics under the direction of an 
editorial board composed of members of the 
Institute. It will be the policy of the Journal to 
print contributions of scientific and technical 
interest to members. In general, it will provide 
an opportunity for members to publish articles 
disclosing new knowledge and new applications 
in the field of the aeronautical sciences. The 
Journal will print a few complete professional 
papers in each issue and a larger number of 
shorter notes and critical comments. The Journal 
will restrict itself to the aeronautical aspect of 
advances in science and engineering. New knowl- 
edge in physics, metallurgy, meteorology, etc., 
having aeronautical applications will be reported 
and discussed in the Journal. Papers read at 
meetings of the Institute will be presented in 
full or digested for the benefit of members who 
are unable to attend. It is hoped to include 
reviews of aeronautical literature and, as soon 
as practicable, a digest of current aircraft books 
and magazines will form a section of the publica- 
tion. The news of the Institute will be furnished 
members through the Journal. The Journal will 
be published quarterly. 

From time to time, the Institute will, when 
it appears desirable, provide members with re- 
prints of articles or other publications of interest. 

An Annual is published giving the membership 
list with the constitution and by-laws of the 
Institute. 


MEETINGS 


Meetings of the Institute are held as often as 
interesting speakers or important aeronautical 
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papers are available. At these meetings, discus- 
sion of the subjects presented is invited. 

Distinguished foreign aeronautical visitors are 
received by members at the club rooms at in- 
formal gatherings. In this way an opportunity is 
provided to reciprocate courtesies received by 
our members when abroad. 

An Annual Meeting is held in January of each 
year at which papers of scientific importance are 
presented and discussed. A review of progress in 
the aeronautical sciences during the year is also 
presented so that members may be kept abreast 
of current activities. 


MEMBERSHIP 
The membership of the Institute consists of: 


Scientific Members 
Engineering Members 
Pilot Members 
Industrial Members 
Technical Members 


Honorary Fellows 
Honorary Members 
Benefactors 
Fellows 


The following qualifications are used in considering the 
applications for memberships: 


Honorary Fellows of the Institute shall be elected by 
the Council and shall be persons of pre-eminence in aero- 
nautics. 


Honorary Members of the Institute may be elected by 
the Council for the period of one year. 


Benefactors shall be those persons who are interested 
and active in the advancement of aeronautics and who 
shall have contributed to the Institute, at one time, the 
sum of not less than one thousand dollars. 


Fellows are elected by the Council from among Members 
of the Institute of at least six months standing. A Fellow 
shall have attained a position of distinction and shall have 
made notable and valuable contributions in one of the 
aeronautical sciences or aeronautical engineering. 


Scientific Members shall be persons of recognized stand- 
ing in one of the aeronautical sciences. 

“Recognized standing”’ when applied to Scientific Mem- 
bers means that applicants should ordinarily be over thirty 
years of age. They should have been engaged for a sub- 
stantial period of time in research, teaching or the applica- 
tion of science to aeronautics. Preferably, this work should 
not be of an immediate commercial application. Recog- 
nized standing is indicated by publications, membership 
in scientific societies, academic position, record of contribu- 
tions to knowledge or the like. 
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Engineering Members shall be persons of recognized 
standing in aeronautical engineering. 

“Recognized standing’’ when applied to Engineering 
Members means that applicants should ordinarily be over 
thirty years of age and, as aeronautical engineers, under 
general direction or independently, have taken responsi- 
bility for important aeronautical engineering work. They 
should have been in active practice of their profession for 
at least five years. Alternatively, they should be persons 
who have been regularly engaged in aeronautical work or 
closely allied work for at least five years, who have, by 
invention or proficiency in aeronautical science, the aero- 
nautical arts, or aeronautical literature, attained a standing 
equivalent to that required of the aeronautical engineers 
as above described. Standing is indicated by publications, 
membership in engineering societies, record of aeronautical 
designs created, and the like. 


Pilot Members shall have qualified as pilots and have 
acquired a recognized standing in the aeronautical field. 

“Recognized standing’’ when applied to Pilot Members 
means pilots whose work has been primarily related to 
flight in its scientific or engineering aspects such as test 
flying. This grade may also include persons of distinction, 
who are pilots but who would not properly be eligible as 
Industrial, Engineering or Scientific Members. 


Industrial Members shall have been engaged in the 
application or development of aeronautics and shall have 
acquired a recognized standing in the aeronautical field. 
Persons whose work is principally executive or commercial 
are classed as Industrial Members. Their standing is 
recognized by their position in the industry. 


Technical Members shall be persons who wish to affiliate 
with the Institute and are engaged in aeronautical engi- 
neering, research, teaching, or study, or scientific or engi- 
neering work in one of the sciences allied to aeronautics. 
Pilots and persons engaged in the aeronautical industry 
or related industries who are interested in the advancement 
of the aeronautical sciences may be eligible as Technical 
Members. Persons who may become eligible for other 
grades of membership upon receipt of more complete infor- 
mation regarding special qualifications may be initially 
classed as Technical Members. 

Graduation from an aeronautical engineering course of 
a school of recognized standing shall be considered the 
equivalent of two years aeronautical experience. 


Foreign Members—Members, of any grade, not residing 
in the United States are considered to be foreign corre- 
spondents with whom the Institute desires to be associated. 

Only Fellows, Scientific Members and Engineering 
Members shall be entitled to vote by letter ballot or at 


meetings of the Institute, but Honorary Fellows, Bene- 
factors, Pilot Members and Industrial Members also may 
serve as officers and on committees. 


Abbreviations—The approved abbreviations of the title 
of the various grades of membership shall be as follows: 


Honorary Fellow............. Hon. F. I. Ae. S. 
Scientific Members...............M. I. Ae. S. 
Engineering Members............. M. I. Ae. S. 


Pilot Members, Industrial Members and Technical 
Members shall not be entitled to the use of the initials 
after their names. 


The entrance fee for all members except Honorary 
Fellows, Honorary Members, Benefactors and Technical 
Members, shall be five dollars and that for Technical 
Members three dollars. 


The annual dues shall be: 


None 
Honorary Members..................... None 
Scientific Members....................-$5.00 
Engineering Members................... $5.00 
Industrial Members... ...... $5.00 
Technical Members. $3.00 


APPLICATION FOR MEMBERSHIP 


Application forms for membership may be 
secured from any member of the Institute or by 
requesting them from the Secretary, 5431 RCA 
Building, Rockefeller Center, New York. 

Applicants are required to furnish the names 
of several persons, preferably members, who are 
acquainted with their professional work. After 
these references have been communicated with 
by the Institute, the Admissions Committee 
recommends the grade of membership. This 
recommendation is then passed on by vote of 
the Council after which an invitation to join 
may be extended to the applicant. 
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The Skyport 


THE R. C. A. BUILDING 
ROCKEFELLER CENTER 
NEW YORK 

THE Skyport, the clubroom and headquarters 
of the Institute, is located on the fifty-fourth 
floor of the largest office building in the world. 
Arrow indicates the location of THE SKYPORT 
at Rockefeller Center. 


The Skyport, the clubroom of the Institute of 
the Aeronautical Sciences, has been provided for 
the use of members by generous donors. An 
anonymous donor has given the rent for the first 
year and Mr. Grover Loening, a member of the 
Council, has provided the furnishings. Others 
have presented books, models, and equipment. 


The surroundings of the Skyport may be better 
understood from the following description of 
Rockefeller Center: 

Rockefeller Center is a permanent century of 
progress in art, industry and entertainment. 
Probably the greatest concentration of men and 
women ever gathered within a single private 
building development will people Rockefeller 
Center. When these buildings have been com- 
pletely tenanted, the total daily population wi!l 
be more than 215,000. From the RCA Building 
you can see for a radius of fifty miles. You can 
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look north to Westchester, east to Long Island, 
south to the Atlantic and west to New Jersey. 

Two international buildings, the British Em- 
pire Building and La Maison Francaise, are now 
complete, and two others, the Palazzo d'Italia 
and a twin structure, are expected to be shortly 
in process of construction. 

The art program of Rockefeller Center has 
scope and variety of unusual proportions. 
Murals, sculpture and other forms of decoration, 
executed by leading artists, may be found on the 
exteriors and interiors of all the buildings now 
completed. Similar ornamentation will be applied 
to future buildings. 

At night the multicolored lights of the city 
and its suburbs offer a picturesque spectacle. 
The neon lights of Times Square, the strings of 
lights in Central Park, the illuminated windows 
of skyscrapers and the streams of moving head- 
lights on highways and trestles provide an un- 
forgettable panorama. 


Necrology 


August Goldsmith, an Industrial Member of 
the I. Ae. S., died in New York City on May 24, 
1933. He was born in Birkenfeld, Germany, on 
May 28, 1860, and came to the United States in 
1878. He founded the B. G. Corporation, manu- 
facturers of spark plugs, in 1917 and was its 
president until his death. 

Edward Tracy Birdsall, a Scientific Member 
of the I. Ae. S., died at his home in New Rochelle, 
New York, on November 6, 1933, at the age of 
seventy years. He was one of the founders of the 
Society of Automotive Engineers and was its 
first Secretary and Treasurer. He was also a 
charter member of the American Institute of 
Electrical Engineers. He graduated as a Mechan- 
ical Engineer from the Stevens Institute of Tech- 
nology in 1886. He was a designer of internal 
combustion engines since 1900, being one of the 
pioneers in the automotive field. For several 
years he had been interested in the development 
of a Diesel engine for aviation purposes. At the 
time of his death he was consulting engineer of 
the Lawrance Engineering and Research Cor- 
poration. 


Annual Meeting to be held at Columbia Uni- 
versity 


An all day meeting of the Institute will be held 
on Wednesday, January 31, 1934, at the Physics 
Building, Columbia University, Broadway at 
120th Street, New York City. 

In the morning, commencing at 10 o'clock, the 
application to aeronautics of Meteorology, Ra- 
dio, Metallurgy, and Physics will be discussed. 
A luncheon will be held at 1 P.M. at the Faculty 
Club of Columbia University for all members and 
guests. During the afternoon the subjects for 
discussion will be Aerodynamics, Fuel, Power 
Plants and Instruments. At five o’clock the An- 
nual Meeting of the Institute will be held. Re- 
ports will be made by the officers. New members 
of the Council will be elected and plans for the 
year discussed. 

In the evening a dinner will be held at the 
Faculty Club. Some of the earliest workers in 
aerodynamics, piloting, ballooning, engines and 
other pioneer fields have been invited to give 
their earliest impressions of what they were try- 
ing to do and to recall early memories. 

The Institute was fortunate, last year, in re- 
ceiving an invitation from Columbia University 
to hold its Founders’ Meeting there. The facilities 
of the Physics Building and the Faculty Club 
are so suitable for such a gathering that the 
Council was pleased to accept the hospitality of 
the Physics Department again. The Council of 
the Institute has passed resolutions appreciating 
the courtesy of Professor George B. Pegram, a 
scientific member of the Institute and Professor 
of Physics at Columbia, who has extended the 
invitation for Columbia and acted as host at our 
meetings. 

It was decided by the Council that, due to the 
limited time available for the meeting, the pres- 
entation and discussion of formal scientific 
papers would not be advisable. It was thought 
preferable to have, instead, several of the leading 
specialists in the fields indicated above to discuss 
and answer questions on recent knowledge and 
developments in their particular fields. In this 
way it is hoped to develop an appraisal of the 
aeronautical sciences in general. With the assis- 
tance of these specialists, the sessions will be 
conducted by President Hunsaker. 
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Some of the designs submitted in the competition for a Seal for the Institute. The award was made to a combination 
of the two seals in the upper left hand corner. See front cover. 
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The Institute Seal 


The Institute needed a seal and the Council 
authorized a competition for designs which might 
be suitable. Forty-six different sketches were sub- 
mitted. There were so many which were excellent 
that it was difficult to make a choice. Two of the 
designs seemed to meet the requirements of the 
competition more nearly than any others. All 
competitors were given the following instruc- 
tions: 

“Tt would be well if aeronautical sciences, as 
distinguished from engineering, could be empha- 
sized. Simplicity and dignity are desired.” 

The final award for the seal which appears on 
the cover of this issue was made to Mr. George 
W. Lewis, Director of Aeronautical Research of 
the National Advisory Committee for Aero- 
nautics and Mr. W. S. Simpson of the Bureau 
of Aeronautics, Navy Department. By combin- 


ing features of both designs, the Council believes 
that it has secured a seal which is representative 
of the underlying purpose of the Institute. The 
Council is appreciative, not only of the designs 
submitted by the winners of the award, but to 
all those who submitted other designs. So that 
members of the Institute can see the various 
sketches, some are reproduced in this issue. 

The Council wishes to thank the following for 
taking part in this competition. Robert A. 
Kaiser, B. Kutakoff, M. P. Weihmiller, H. E. 
Weihmiller, V. Schreckengost, D. J. Leisk, H. 
McCallister, H. K. Growald, D. E. Dunlap, W. 
Klemperer, K. M. Daland, G. F. Cummings, H. 
Case, J. A. Roche, S. T. Payne, J. L. Fechter, 
R. I. Markey, R. Funk, M. Cozzoli, E. L. 
Westdal, P. Pevney, J. R. Douglass, K. Huerttle. 


Donots to The Skyport 


The Institute has received gifts from donors 
to whom it wishes to express its grateful appre- 
ciation. 

To the anonymous donor who has provided 
headquarters and a clubroom at Rockefeller 
Center for one year the Institute is under great 
obligation. To the owners of the great building 
development, Rockefeller Center, Inc., the Insti- 
tute is indebted for many favors. It is hoped that 
Rockefeller Center will become the center of 
aeronautical activity in New York. The various 
Curtiss-Wright-Sperry corporations and _ the 
Manufacturers Aircraft Association have located 
in the RCA Building. 

On October 1st, when the Institute signed its 
lease it did not own any furnishings. Mr. Grover 
Loening volunteered to furnish and decorate the 
office and clubrooms. By this donation, which 
amounted to over two thousand dollars, Mr. 
Loening becomes the first Benefactor of the 
Institute, according to the requirements of the 
By-Laws. 


Other donors to whom the Institute wishes to 
express its thanks for gifts are as follows: 

Aluminum Company of America—Aluminum 
chairs. 

Mr. Luis de Florez—Radio receiver. 

Major Lester D. Gardner—Aeronautical li- 
brary of several hundred volumes, file of aero- 
nautical pamphlets and photographs. 

Dr. Henry A. Gardner—Laboratory manual 
“Paints, Varnishes, Lacquer, and Colors.” 

Mr. Edward J. Steichen—Photographic mural 
for the Skyport. 

R.K.O. Corporation—Permission to use Mr. 
Steichen’s mural which was originally made for 
the Center Theater, New York. 

Prof. Alison Reppy—Three volumes of ‘‘Air 
Law Review.” 

Mr. T. P. Wright—Two Curtiss airplane 
models. 

Mr. R. H. Fleet—Typewriter and airplane 
model. 

Mr. Paul Litchfield—Photographic Mural of 
airship ‘‘Macon.” 
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The photographic mural by Edward Steichen showing development of flight. 
Models on exhibition during week of opening. 


Aluminum is used in THE Sxyport for furnishings and decorations. 
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oo Members’ writing and reading room at THE SKYPORT. 


